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Abstract 

We present a review of gravitating particle-like and black hole solutions with non- 
Abelian gauge fields. The emphasis is given to the description of the structure of 
the solutions and to the connection with the results of flat space soliton physics. 
We describe the Bartnik-McKinnon solitons and the non-Abelian black holes arising 
in the Einstein- Yang-Mills theory, and consider their various generalizations. These 
include axially symmetric and slowly rotating configurations, solutions with higher 
gauge groups, A-term, dilaton, and higher curvature corrections. The stability issue 
is discussed as well. We also describe the gravitating generalizations for fiat space 
monopoles, sphalerons, and Skyrmions. 
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1 Introduction 



Non-Abelian solitons play an important role in gauge theories of elementary 
particle physics [94,3,267]. In most cases, the effects of gravity can safely be 
neglected, perhaps except for the very heavy magnetic monopoles in some 
Grand Unified Theories. Besides, it has recently become clear that solitons 
are equally important also in string theory, where gravity is essential. A large 
number of solutions for such gravitating lumps have been obtained [118], how- 
ever, almost all of them are Abelian. At the same time, the gauge group in 
string theory is quite complicated, which suggests that non- Abelian solutions 
can be important as well. From the point of view of General Relativity (GR), 
the theory of gravitating non- Abelian gauge fields can be regarded as the most 
natural generalization of the Einstein-Maxwell (EM) theory. It is therefore rea- 
sonable to study gravitating gauge fields and, in particular, to check whether 
the standard electrovacuum results have natural generalizations. In view of 
this, the aim of the present review is to discuss non- Abelian gravity-coupled 
solitons and black holes. The concepts of solitons and lumps refer in this text 
to any particle-like solutions of a non-linear field theory. Such solutions are 
asymptotically flat, topologically trivial and globally stationary, although not 
necessarily stable. 

The first example of gravitating non- Abelian solitons was discovered by Bart- 
nik and McKinnon (BK) in the four- dimensional Einstein- Yang-Mills (EYM) 
theory for the gauge group SU(2) [19]. This example can be regarded as canoni- 
cal in the sense that solutions in other models of gravitating non- Abelian gauge 
fields studied so far inevitably share a number of characteristic common fea- 
tures with the BK particles. Soon after the BK discovery it was realized that, 
apart from solitons, the EYM model contains also non-Abelian black holes 
[324,215,33]. As these manifestly violate the non-hair conjecture, they have 
attracted much attention and stimulated a broad search for other black hole 
solutions in models of four-dimensional Einstein gravity with non-linear field 
sources. The results obtained have led to certain revisions of some of the basic 
concepts of black hole physics based on the uniqueness and no-hair theorems. 
Specifically, it has been found that the violation of the no-hair conjecture is 
typical for gravitating non-Abelian gauge theories, especially for those models 
which admit solitons in the flat spacetime limit. 

Let us recall (see [164] for a more detailed account) that the existing classi- 
fication of black holes in GR is based on Hawking's strong rigidity theorem, 
stating that a stationary black hole is either static (with non-rotating horizon) 
or axisymmetric [159]. This theorem uses only fairly general assumptions like 
the weak energy condition (although some of them are being critically revised 
[91,90]). For vacuum and electrovacuum black holes with non- vanishing surface 
gravity at the horizon, Israel's theorems [176,177] ensure that staticity implies 
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spherical symmetry. It follows then that the static Schwarzschild and Reissner- 
Nordstrom (RN) solutions, respectively, are unique. For zero surface gravity, 
static electrovacuum black holes belong to the Majumdar-Papapetrou family 
[166]. In the stationary chain of uniqueness theorems asserts that the 

regularity of the event horizon and the asymptotic behaviour specified by mass, 
angular momentum and Coulomb charges determine the solutions completely. 
As a result, stationary electrovacuum black holes are necessarily axisymmet- 
ric and should belong to the Kerr-Newman family [77,243,75]. Apart from the 
vacuum and electrovacuum cases, uniqueness has been established for super- 
gravity models containing in four- dimensions sets of scalars parameterizing 
coset spaces and the corresponding multiplets of the U(l) fields [62,60]. 

Another well-known statement, the no-hair conjecture, claims that the only 
allowed characteristics of stationary black holes are those associated with the 
Gauss law, such as mass, angular momentum and electric (magnetic) charges 
[277,21]. Thus, for example, it follows that black holes cannot support ex- 
ternal scalar fields, since there is no Gauss law for scalars. On the other 
hand, black holes with a spin 3/2- field, say, are allowed, since they carry a 
conserved fermion charge given by the Gauss flux integral [6]. The conjec- 
ture has been proven for non-interacting boson [25,24,23] and fermion fields 
[157,303,304] of various spins, as well as for some special non-linear matter 
models [4,22,242,162,163,298]. 

Although the uniqueness and no-hair theorems had actually been proven only 
for special types of matter, the appealing simplicity of these assertions inspired 
a widespread belief in their possible general validity. However, after the discov- 
ery of the EYM black holes it became clear that this is not the case. First, the 
EYM black holes possess a short-ranged external non-Abelian gauge field and 
are not uniquely specified by their mass, angular momentum and conserved 
charges. The no-hair conjecture is therefore violated. Second, static EYM black 
holes with non-degenerate horizon turn out to be not necessarily spherically 
symmetric [200,202]. This shows that Israel's theorem does not generalize to 
the non-Abelian case. A similar phenomenon has also been observed in the 
EYM-Higgs model [273]. Next, the perturbative considerations suggest that 
non-static EYM black holes with rotating horizon do not necessarily have 
non-zero angular momentum [327,66], which shows that the Abelian static- 
ity conjecture [78] does not straightforwardly apply either [170,300]. In addi- 
tion, the Frobenius integrability conditions for two commuting Killing vectors 
are not automatically fulfilled for self-gravitating Yang-Mills fields [170,165], 
and so the Ricci circularity condition is not guaranteed. The standard Lewis- 
Papapetrou parameterization of a stationary and axisymmetric metric used 
for the uniqueness arguments can therefore be too narrow. However, explicit 
examples of the circularity violation are not known yet. Finally, the three- 
dimensional reduction of the Yang-Mills action in the presence of a Killing 
symmetry does not lead to the standard sigma-model structure essential for 
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the uniqueness proof [161,126]. 



All this shows that a number of very important features of electrovacuum black 
hole physics cease to exists in the EYM theory, and this seems to happen gener- 
ically for models with non-Abelian gauge fields. At the same time, what might 
seem surprising from the traditional point of view, finds a natural explanation 
in the context of flat space soliton physics [94,3,267]. For example, the exis- 
tence of hairy black holes can often be directly inferred from the existence of 
solitons in flat spacetime. For small values of Newton's constant the implicit 
function theorem ensures that flat space solitons admit weakly gravitating, 
globally regular generalizations - weak gravity can be treated perturbatively. 
It turns out that for a large class of matter models the same argument can be 
used to show the existence of solutions in a small neighbourhood of the point 
r h = in parameter space, where r h is the event horizon radius [188]. This 
implies that a weakly gravitating lump can be further generalized to replace 
the regular origin by a small black hole. The result is a hairy black hole whose 
radius is considerably smaller than the size of the lump surrounding it. One 
can argue in the spirit of the no-hair conjecture that a small black hole can- 
not swallow up a soliton which is larger than the black hole itself. However, 
a larger black hole must be able to do this, and indeed, the radius of a hairy 
black hole inside a soliton usually cannot exceed some maximal value. 

In view of the arguments above, the most surprising fact is the existence of the 
BK particles, which have no flat space counterparts, and of their black hole 
analogues, whose radius can be arbitrary. The pure EYM theory is character- 
ized by the fact that all fields are massless, and the solutions of the theory 
can therefore be thought of as bound states of two non-linear massless fields. 
As a result, the EYM solitons and black holes exhibit a number of special fea- 
tures which distinguish them from solutions of other gravitating non-Abelian 
gauge models. For example, the BK particles provide the only known example 
of solitons, whether gravitating or not, which admit regular, slowly rotating 
generalizations [69]. Solitons of other models, such as the t'Hooft-Polyakov 
monopoles, say, do not rotate perturbatively [66,171], which means that they 
either cannot rotate at all or their angular momentum assumes discrete values. 
The EYM black holes are uniquely distinguished by their peculiar oscillatory 
behaviour in the interior region [115,59], which is reminiscent of some cosmo- 
logical models [27]. 

On the other hand, some features of the EYM solutions are shared by those of 
other non-linear models. The BK solitons, for example, exhibit a remarkable 
similarity [134] with the well-known sphaleron solution of the Weinberg-Salam 
model [208], which allows one to call the BK particles EYM sphalerons. (Let 
us remind to the reader that the term "sphaleron" refers to the static saddle 
point solution in a gauge field theory with vacuum periodicity [179]. Sphale- 
rons are characterized by half-integer values of the Chern-Simons number of 
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the gauge field and can be reduced to any of the nearest topological vacua 
via a sequence of smooth deformations preserving the boundary conditions 
[237,208]. Such objects are likely to be responsible for the fermion number 
non-conservation at high energies and/or temperatures [217]). Other prop- 
erties of EYM solitons and black holes, such as their nodal structure and 
discrete mass spectrum, are generic and shared by practically all known solu- 
tions in models with gravitating non-Abelian gauge fields. Such solutions can 
be thought of as eigenstates of non-linear eigenvalue problems, which accounts 
for the discreteness of some of their parameters. 

The list of the non-linear models with gravitating non-Abelian gauge fields, 
investigated during recent years, contains, apart form the pure EYM the- 
ory, also various generalizations. These include the dilaton, higher curvature 
corrections, Higgs fields, and a cosmological constant. The Einstein- Skyrme 
model has also been studied in detail. All these theories admit gravitating 
solitons and hairy black holes. The pure EYM theory has been studied most 
of all, some of its features being typical for all other models. For this reason, 
the central part of our review will be devoted to a description of the basic 
EYM solutions as well as their direct generalizations. However, other impor- 
tant solutions, such as the gravitating monopoles and Skyrmions, will also be 
duly described. We will concentrate only on the most important results, but 
our list of references is quite complete. Several review articles on the related 
subjects are available [140,181,254,37,232,165,21,161], but overlaps with the 
present text are small. 

The plan of the paper is as follows. In Sec. 2 we give the basic definitions, 
briefly discuss symmetries of the non-Abelian gauge fields, and derive the 
reduced two-dimensional Lagrangians for the spherically symmetric models. 
Sec. 3 is devoted to the BK solutions, their sphaleron interpretation and various 
generalizations. Basic properties of the EYM black holes are discussed in Sec. 4. 
Sec. 5 and Sec. 6 contain a discussion of stability of EYM solitons and black 
holes and an analysis of their rotational excitations. Gravitating solutions in 
models admitting solitons in flat space are discussed in Sec. 7. These include 
the t'Hooft-Polyakov monopoles, the Yang-Mills-Higgs sphalerons, and the 
Skyrmions. Some other important results not mentioned in the main text are 
briefly discussed in Sec. 8. Most of the solutions described below are known only 
numerically, although their existence has been established rigorously in some 
cases. For this reason we present a number of figures and tables describing 
the numerical results. All of them have been produced with the "shooting to 
a fitting point" numerical procedure, described in [266]. 
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2 General formalism 



In this chapter the field-theoretical models, whose solutions will be consid- 
ered in the next chapters, are introduced. These are the EYM theory and its 
generalizations, including the dilaton, the Higgs fields, as well as the Einstein- 
Skyrme model. We briefly discuss symmetry conditions for gauge fields and 
derive the effective two-dimensional Lagrangians for all models in the spheri- 
cally symmetric case. 



2.1 Einstein- Yang- Mills theory 

The basic model which will be discussed below is the four-dimensional EYM 
system for a compact, semi-simple gauge group Q. The Lie algebra of Q 
is characterized by the commutation relations [T , T 6 ] = if a bcXc (°, b,c — 
1, . . . , dim(£)). The basic elements of the model are (M,g^ u , A), where M. is 
the spacetime manifold with metric g^, and the Lie-algebra- valued one- form 
is A = A^dx^ = T a A®dx^. We choose the standard action 

s ™ = /{"Tab* " W^ tlF - F ""} (2 ' 1) 

where g is the gauge coupling constant, = d^A u — d^A^ — iL4 M , A v ] = 
T a F£ v , and K > is a normalization factor: tr (T a T 6 ) = K 5 ab . For £=SU(2) 
we choose T a = |r a , K = 1/2, f abc = e abc with r a being the Pauli matrices. 

We adopt the metric signature (H ), the spacetime covariant derivative 

will be denoted by V, and the Riemann and Ricci tensors are R 01 ^ = — 
. . . and = R a ^ au , respectively. The gravitational constant G is the only 
dimensionful quantity in the action (the units h — c — 1 are understood). 

Apart from the general spacetime diffeomorphisms, the EYM action is invari- 
ant with respect to the gauge transformations of the gauge field A: 

A^-UOV + i^lT 1 , F^-U^IT 1 , g,u^9^, (2.2) 

where U(x) G Q. In addition, the Yang-Mills (YM) part of the action displays 
the conformal symmetry 

g^ -> Sl(x) g^ , A^ -> Ap, F^ -> F^, (2.3) 
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which is not, however, a symmetry of the Einstein-Hilbert action. The varia- 
tion of the action (2.1) gives the Einstein equations 

R^-^Rg^ = 8nGT^ (2.4) 

with the gauge-invariant YM stress-energy tensor 

T, v = tr (-F^FS + \ 9ilv F a ?F a ^ , (2.5) 

and the YM equations 

= , (2.6) 

where the gauge covariant derivative is = V M — ]. Owing to the 

conformal invariance the stress tensor is traceless, = 0. As a consequence of 
the gauge invariance, the dual field tensor *F tlv = (l/2)y/^g e ilvpa F ( " J (e 0123 = 
1) satisfies the Bianchi identities 

£> M *F^ = 0. (2.7) 

Asymptotically flat solutions in the theory have well-defined ADM mass and 
angular momentum [11,1]. Similarly, one can define conserved Lie-algebra- 
valued electric and magnetic charges [339,2]. Note that the straightforward 
definition [339] 




where the integration is over a two-sphere at spatial infinity, is not gauge in- 
variant in the non-Abelian case [281,92], and hence requires the gauge fixing. 
This reminds of a similar situation for the ADM mass, which was originally 
defined only for distinguished coordinate systems [11]. On the other hand, in 
the presence of a timelike Killing symmetry the ADM mass can be covariantly 
expressed by the Komar formula [209] . For the YM charges there exists a sim- 
ilar gauge invariant construction [227] using the Lie-algebra valued "Killing 
scalar" K, where D^K = 0. Now, it is important that these invariant defini- 
tions can be generalized to the case where the configuration has no symmetries 
at all. Indeed, the symmetries always exist in the asymptotic region, and these 
can be used in order to put the surface integrals into the covariant form. As a 
result one can define the ADM mass [1] and conserved gauge charges [2] for an 
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arbitrary isolated system in the completely covariant and gauge- invariant fash- 
ion. For example, for the gauge field given by A = A + a, where A = ill dU , 
and a — > as r — > oo, one can define = D^ay — -Dj,a M . Here D is the 
covariant derivative with respect to A. The conserved and gauge invariant (up 
to global gauge rotations) charges Q a and P a are then given by (2.8) with F 
replaced by tr^FTjT 1 ) [2]. 

2.2 Spacetime symmetries of gauge fields 

The issue of symmetry for non-Abelian gauge fields has been extensively stud- 
ied in [274,29,125,154,156,155,236,150,160,99,20,211,210,17,213,70,72,18]. The 
symmetry conditions for gauge fields in the infinitesimal language were first 
formulated in [29,125]. An EYM field configuration will respect spacetime 
isometries generated by a set of Killing vectors £ m if the metric is invariant, 

A m ^ = 0, (2.9) 

under the action of the isometry group, while the corresponding change in the 
gauge field, A^ — > A^ — e m C^ m A^, can be compensated by a suitable gauge 
transformation, A^ — > A^ + e m D^W m , such that [29,125] 

C u A, = D^W m . (2.10) 

Here C stands for the Lie derivative. In order to actually find A^ for a given set 
of Killing vectors £ m the procedure is first to solve the integrability conditions 
for (2.10), which gives the WJs, and then to solve (2.10) for A^ [125,20]. 

Throughout this article we will assume the 2 + 2 block-diagonal form of the 
metric 

ds 2 = g a pdx a dx p + h MN dx M dx N , (2.11) 

where x a = {a; ,^ 1 }, x M = {x 2 ,x 3 }, and g a p depend only on x a . In the case 
of spherical symmetry the Killing vectors £ m generate the so (3) algebra, and 
the most general solution to (2.9) can be parameterized as 

ds 2 = g a(3 dx a dxP - R 2 (dd 2 + sin 2 &d<p 2 ), (2.12) 

where R depends on x a = {t,r}. The corresponding most general solution 
to Eq.(2.10) for the gauge group SU(2) is sometimes called Witten's ansatz 
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[103,335,125,17]: 

A = a T r + i(l - Re iu) [T r , dT r ] +lmw dT r . 



(2.13) 



Here the real one-form a = a a dx a = a c?t + a r dr and the complex scalar w 
depend only on x a . The position-dependent gauge- group generators 

T r = n a T a , T^ = ^T r , T v = -^—d^T r , (2.14) 

sm v 

where n a = (sin -d cos ip, sin-^siny?, cos-#), obey the standard commutation re- 
lations, such that [T r ,T$] = iT^. The ansatz (2.13) is invariant under the 
U(l) gauge transformations generated by U= exp(i f3(t, r) T r ), under which 

a a -> a a + d a/ 3, u> -> e l/ V (2.15) 

i.e., a Q and w transform as a two-dimensional Abelian vector and a com- 
plex scalar, respectively. The four independent real amplitudes in (2.13) con- 
tain therefore one pure gauge degree of freedom. One can choose and 
(arg w) as the three gauge-invariant combinations. Under parity 
transformation, -t? — > 7r — p ^ n + <p, one has T r — > — T r , T,g — > T#, 
— > — T v . As a result, the effect of parity on (2.13) is 

P: a a -> -a a , w ^ (2.16) 

and hence among the four real amplitudes one is parity-even, and three are 
parity-odd. 

After the gauge transformation with U=exp(i T 2 1?) exp(i T 3 </?) [9], an equiv- 
alent form of the ansatz (2.13) in terms of the constant group generators T a 
can be obtained: 

A = a T 3 + Im (wT+) dd - Re (wT+) sin $ dp + T 3 cos dp, (2. 17) 

where T + = Ti + iT 2 . The residual gauge freedom (2.15) is now generated 
by U= exp(i /3(t, r) T 3 ). The expression (2.17) is sometimes easier to use than 
(2.13), but it contains the Dirac string singularity. Another frequently used 
gauge is w = w*, in which case the gauge field admits the following useful 
representation: 

1 — vo 

A = aT r + i UdlT 1 , where U = exp(i7rT r ). (2.18) 

2 

The corresponding field strength reads 
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F = - T r da - i dw A [T r , dT r ] 

- w a A dT r + T r (w 2 - 1) d-d A sin dp. 



(2.19) 



Note that the components F aip and -F Q do not vanish. 

Unlike the situation in the Abelian case, in order to achieve the temporal 
gauge condition ao = one has to make a gauge transformation which renders 
the whole configuration time- dependent. Therefore, for static fields one cannot 
set a = 0, unless f2 = a — (arg w) - =0. The second gauge- invariant com- 
bination, Q r = a r — (arg w)', vanishes in the static case by virtue of the YM 
equations. As a result, one can always reduce a r to zero by a gauge transfor- 
mation, and hence the most general static, spherically symmetric SU(2) YM 
fields can be parameterized by two real functions: ao and itQ. 

In the static, purely magnetic case the SU(2) ansatz (2.17) with a = Imw = 
can be generalized to the gauge group SU(N) [333,213,18]. Such a generaliza- 
tion includes instead of one function w, N — 1 independent real amplitudes 
Wj(r), j — 1 . . . N — 1. The gauge field potential A is given by a Hermitean 
N x N matrix, whose non-vanishing matrix element are 

Ajj = - (N - 2j + 1) cos dtp, j — 1...N, 
A j>j+1 = (A j+1>j )* = l Wj e, j = 1 . . . N - 1. (2.20) 

Here = y~j(N — j) (i d$ + sin dip) and the asterisk denotes complex con- 
jugation. 

In the static, axially symmetric case, the metric is given byQ 

d s 2 = e s dt 2 - e c dp 2 - e^(dp 2 + dz 2 ), (2.21) 

where 5, (, and \i depend on x a = {p, z}. The purely magnetic SU(2) gauge 
field reads [269] 

A = a + {Re w T p + (Im w - v) T 3 }d<p , (2.22) 



In systems with a Higgs field the YM equations do not always imply that a r 
vanishes in the static case. The most general spherically symmetric YM fields include 
then ao, a r , and w = w* . 

2 This metric fulfills the circularity condition and gives rise to non-trivial EYM 
solutions. It is unclear, however, whether other, non-circular solutions can exist as 
well. 
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where the one-form a = a p dp + a z dz and the complex scalar w depend on p 
and z. The group generators are given by 

T p = cosz/^T! + sinz/^T 2 , T^ = v~ 1 d v T p , (2.23) 

with integer v. The field (2.22) also has the residual U(l) invariance expressed 
in the form (2.15) and generated by U= exp(i /3(p, z) T^). 

Note that the ansatz in Eq.(2.17) can be generalized to describe spherically 
symmetric YM fields with arbitrary winding numbers: 

A = a T 3 + Im (wT+) dti-vRe (wT + ) sin &d(p + vT 3 cos dip. (2.24) 

Here an integer v has the meaning of the Chern number of the U(l) bundle 
[17]. However, for v ^ 1 the YM equations impose the following condition for 
non- vacuum fields: w — 0. This implies that all configurations with v ^ 1 are 
embedded Abelian. Since we are interested in non-Abelian solutions, for which 
w does not vanish identically, we shall always assume that v — 1 in (2. 24)0. 

2.3 Dimensional reduction and scalar fields 

We shall mainly be considering below spherically symmetric systems. It is 
convenient to derive the field equations in the spherically symmetric case by 
varying the effective two-dimensional action (the Hamiltonian approach to 
the problem was discussed in [100,28]). Let us insert into the four-dimensional 
(4D) EYM action (2.1) the ansatz (2.13) for the gauge field and the metric 
parameterized as 

ds 2 = a 2 N(dt + adrf - — dr 2 - R 2 (d$ 2 + sin 2 $tp 2 ), (2.25) 

with a, N, a, and R depending only on x a = {t,r}. Integrating over the 
angles and dropping the total derivative, the result is the two-dimensional 
(2D) effective action [68] 

S = J {rri - am + L m }ad 2 x, (2.26) 

where m is defined by the relation 

(dR) 2 = d a Rd a R = ^ - 1. (2.27) 

R 

3 The change v — > — v in (2.24) can be achieved by a gauge transformation. 
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From now on the coordinates x a and all other 2D quantities are assumed to 
be dimensionless, after the rescaling x a — > x a /L, where L is a length scale. 
The dimensionless gravitational coupling constant k in (2.27) is proportional 
to G/L 2 . For the scale-invariant YM theory L is arbitrary. It is convenient to 
choose L = y/AnG/g, in which case k — 1. The matter field Lagrangian L m in 
(2.26) then reads 

Lym = ~ Uf a " + \Dw\ 2 - ^ (\w\ 2 - l) 2 , (2.28) 

where f a p = d a ap — dpa a and D a = d a — ia a are the 2D Abelian field strength 
and the gauge covariant derivative, respectively. One can see that the residual 
U(l) gauge symmetry (2.15) arises naturally at the level of the reduced action. 

Note that, since the theory under consideration is purely classical, the EYM 
length scale L and the corresponding mass scale, M—kL/G, do not contain 
Planck's constant. Restoring for a moment in the formulas the speed of light, 
c, one has 

It is worth noting that, up to the replacement \J\~k j 1 g — > e, this corresponds 
to the classical system of units based on G, c, and the charge e that had 
been introduced in [182] before Planck's units were discovered. Dividing and 
multiplying by Vhc, one obtains L = y/a L P \ and M=^ l /a M P i, where Lpi and 
Mpi are Planck's length and Planck's mass, respectively, and a = An/hcg 2 . 
Since one can expect the value of \fa not to be very far from unity, this shows 
that the units (2.29) are actually closely related to Planck's units. 

a) Let us consider generalizations of the EYM theory with additional scalar 
fields. The first one is the EYM-dilaton (EYMD) model with the 4D matter 
action 

Symd = / |^ (V^-^e^tr V-jv^^x, (2.30) 

where 7 is a parameter. In the spherically symmetric case, with <fi = (f>(t, r), the 
4D EYMD action reduces to the 2D form (2.26) with the matter Lagrangian 
L m given by 

iYMD = Y (^) 2 + ^^YM, (2.31) 
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where L YM is specified by (2.28). The length scale is now fixed, L = y/AnG/g, 
and k — 1. 

b) Next, consider the EYM-Higgs (EYMH) model with the Higgs field $ = 
$ a T a in the adjoint representation of SU(2). The 4D matter action is the sum 
of the YM term and 

S u = J tr - ^ ($ a $ a - v 2 ) 2 ^ d 4 x , (2.32) 



where = <9 M $ — $]. The spherically symmetric Higgs field is 

$ = v<pT r , (2.33) 

= (f)(t,r) being a real scalar. This gives the 2D matter Lagrangian in (2.26) 
L m = L YM + ^ («90) 2 - - j R 2 (0 2 - I) 2 - (2.34) 

Here the length scale is L = 1/Mw, e = Mh/a/2Mw, where Mw = gu and 
M H = v^Af are the vector boson mass and the Higgs boson mass, respectively. 
The gravitational coupling constant is k = AnGv 2 . 

c) For the EYMH model with a complex Higgs field $ in the fundamental 
representation of SU(2) the 4D Higgs field action is 

Sh = J ({D^y D"$ - ^ ($t$ - v 2 ) 2 ^ ^d A x , (2.35) 

where -D M $ = <9 M $ — iA^§. The spherically symmetric Higgs field is given by 

$ = vcj) exp(i£T r )|a) (2.36) 



where = (f>(t,r) and £ = r) are real scalars, while |a) is a constant unit 
vector of the 2-dimensional complex vector space, (a\a) = 1. The correspond- 
ing 2D matter field Lagrangian is 

L m = L YM + R 2 \Dh\ 2 S\ 2 -^R 2 (0 2 - l) 2 , (2.37) 



where h = 0exp(i£/2) and D a = d a — ia a /2, the other parameters being the 
same as in the triplet case. 
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d) Finally, we shall consider the SU(2) Skyrme model. The 4D matter action 
is 

Ssk = /tr j-^iV + ^ ( V") J V=9d 4 x, (2.38) 

where A M = Wd^U and = L4 M , A v \. The spherically symmetric chiral field 
is specified by 

U = exp(2^(t,r)T r ), (2.39) 

which leads to the 2D Lagrangian 

Lm=(j+ sin 2 X ) {d X f -{R 2 +\ ^ x) (2.40) 

The length scale is L = 1/ef, the gravitational coupling constant being k = 
AnGf 2 . 

The Einstein-matter coupled field equations are obtained by varying (2.26) 
with respect to N, a, a, and the matter variables in L m . After varying one 
can set the shift function a to zero, but one cannot do this before varying, 
since for time-dependent fields one of the Einstein equations would be lost 
[68]. In the static case, on the other hand, one can get rid of a from the very 
beginning. If the Schwarzschild gauge condition R = r can be imposed, which 
is often the case, the Einstein equations for static fields look especially simple. 
Note that the 4D Einstein equations read in the dimensionless notation 

G^ = 2kT^. (2.41) 

For all matter models described above the matter Lagrangian in the static, 
purely magnetic case reduces to 

L m = -(NK + U), (2.42) 

where K and U depend only on the matter variables, and N = 1 — 2nm/r. 
As a result, varying the reduced action (2.26) with respect to m and a gives 
the independent Einstein equations 

rri = NK + U = r 2 T °, (In a)' = 2kK = k^ (T ° - T r r ), (2.43) 

from which the components of the (dimensionless) stress tensor can be read 
off. The matter field equation are obtained by varying aL m . The dimensionless 
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ADM mass is given by m(oo). The dimensionful mass is (kL/G) m(oo). 

We shall consider below time-dependent problems only for the EYM fields, in 
which case the full system of equations following from (2.26) and (2.28) reads 

d a (r 2 af a ^) = 2alm(wD^w)*, (2.44) 
D a (aD a w) = ^(\w\ 2 -l)w, (2.45) 

(lna)' = l(j^- 2 \Dow\ 2 + \D r w\ 2 ) , (2.46) 
m = 2N Re D w (D r w)*, (2.47) 

m' = -jfapf* + ^lAH 2 + N\D r w\ 2 + ^ (\w\ 2 - l) 2 (2.48) 
with the notation of (2.28) and the asterisk denoting complex conjugation. 
2.4 Static EYM fields. Embedded Abelian solutions 

For static EYM fields the results of this chapter can be summarized as follows. 
The most general spherically symmetric SU(2) gauge field can be parameter- 
ized as 

A = a T 3 dt + w (T 2 dd - Ti sin •& dip) + T 3 cos d dip, (2.49) 

and the spacetime metric is 

ds 2 = a 2 Ndt 2 -^dr 2 - r 2 (dd 2 + sin 2 # dip 2 ). (2.50) 

The amplitudes ao, w = w*, m, N = 1 — 2m/ r, and a depend on r. The EYM 
equations (2.44)-(2.48) assume the form 
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(2.51) 
(2.52) 

/ 2 + ^(^ 2 -l) 2 , (2.53) 

(2.54) 
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In the next two chapters we shall study solutions to these equations. The 
simplest one is 

a = 0, w = ±l, m = M, a = 1, (2.55) 

which describes the Schwarzschild metric and a pure gauge YM field. The 
simplest solution with a non-trivial gauge field describes colored black holes 
[13,87,330,261,184] 

O 2M O 2 + 1 
a = a (oo) + -, w = 0, N = 1 + ^-=— , a = 1. (2.56) 

The metric is RN with the electric charge Q and unit magnetic charge. This 
solution is embedded Abelian in the sense that the SU(2) gauge field potential 
(2.49) is a product of the potential of the U(l) dyon, A = a dt + cosfidip, and 
a constant matrix. The appearance of such a solution is not surprising: since 
the SU(2) gauge group contains U(l) as a subgroup, the set of solutions in 
the EYM theory includes, in particular, the embeddings of all electrovacuum 
solutions. These are often found by integrating the EYM equations. In this 
way, for example, the colored black holes and their generalizations with rota- 
tion and cosmological term [10,189-191], as well as solutions with cylindrical 
and plane symmetries [246,248,247,249,245] have been obtained. Note that 
the embedded Abelian solutions are not completely equivalent to their U(l) 
counterparts. For example, the Dirac string for (2. 49), (2. 56) can be globally 
removed by passing to the regular gauge (2.13) - the SU(2) bundle is trivial. In 
addition, unlike their U(l) counterparts, the colored black holes are unstable 
[221]. 



3 Particle-like solutions 

The rapid progress in the theory of gravitating solitons and hairy black holes 
started after Bartnik and McKinnon had discovered particle-like solutions for 
the EYM equations (2.51)-(2.54) [19]. This came as a big surprise, since it 
is well known that when taken apart, neither pure gravity nor YM theory 
admit particle-like solutions. Indeed, the existence of stationary gravitational 
solitons is ruled out by Lichnerowicz's theorem [224]. One can go farther than 
this and consider regular configurations with finite mass and arbitrary time- 
dependence, as long as they do not radiate their energy to infinity. In par- 
ticular, system which are exactly periodic in time do not radiate, and such 
solutions are also ruled out [139]. It is worth noting, however, that the vacuum 
theory of gravity admits, "quzi-solitons" , geons, which exist for a long time 
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as compared to the characteristic period of the system before radiating away 
their energy [64] Q 

Another famous statement "there are no classical glueballs" asserts that the 
pure YM theory in flat space does not admit finite energy non-singular solu- 
tions which do not radiate energy out to infinity [94,105,260,95,96]. This can 
be readily seen in the static situation, since the spatial components of the 
stress-energy tensor T^ u then satisfy [105] 



This is the consequence of the identity T ik = djiT^Xk) implied by the conser- 
vation law djT- = 0. The physical meaning of (3.1) is that the total stresses in 
an extended object must balance. It follows then that the sum of the principal 
pressures, T*iPi, where pi = —T\ (no summation), cannot have a fixed sign. 
However, for the YM field the scale invariance implies that S^pj = Too > 0, 
which shows that the system is purely repulsive and the force balance is im- 
possible. The argument can be extended to exclude solutions with arbitrary 
time-dependence, as long as they do not radiate [260,95]. It is interesting that 
the geon-type solutions in the pure YM theory are also excluded - any field 
configuration initially confined in some region falls apart in the time it takes 
light to cross the region [96]. This can be viewed as another indication of the 
purely repulsive nature of the gauge field. 

Although for purely attractive or repulsive systems the balance cannot exist, 
the situation changes if the system includes interactions of both types. The 
example is provided by the YM-Higgs models (2.32), (2.35). For the scalar 
Higgs field one has T^pi < 0, which corresponds to pure attraction. However, in 
a combined system with both gauge and Higgs fields repulsion and attraction 
can compensate each other, which leads to the existence of soliton solutions. 
These will be considered in Sec. 7 below. In a similar way, the coupling of the 
repulsive YM field to gravity can lead to a force balance, as is illustrated by 
the BK example. 

It should be stressed, however, that the presence of both attractive and repul- 
sive interactions is necessary but not sufficient for the existence of equilibrium 
configurations. This becomes especially clear when one takes into account a 
number of the no-go results for the EYM system [331,106,234,235,238,121]. 
These, for example, rule out all non-trivial EYM solutions in three spacetime 
dimensions [106], and large classes of solutions in four dimensions. In partic- 
ular, all charged 4D EYM solitons are excluded [121]. In view of this one can 
wonder about the sufficient condition which ensures the existence of the BK 

4 Note though, that no rigorous existence proof for geons has been given so far. 




(3.1) 
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solutions. Unfortunately, this condition is unknown. Even though the exis- 
tence of the BK solutions can be deduced from a complicated analysis of the 
differential equations [293,287,55], it still does not find a clear explanation in 
physical terms. For example, the topological arguments, which proved to be 
very useful in the YMH case, do not apply in the EYM theory. 

It is interesting that, even though the topology is different, there exists a 
very useful analogy between the EYM theory and the flat space doublet YMH 
model. It has already been mentioned that the BK particles are in some ways 
similar to the sphaleron solutions of the Weinberg-Salam model [237]. The 
resemblance is two-fold. First, both sphalerons and BK solitons can be thought 
of as equilibrium states of a pair of physical fields one of which is repulsive (YM 
field) while another one is attractive (gravity or the Higgs field) . Secondly, both 
types of solutions relate to the top of the potential barrier between distinct 
topological vacua of the gauge field. In this sense one can call the BK particles 
EYM sphalerons [134]. 

In view of their importance, we shall describe below in this chapter the BK 
solutions and the corresponding sphaleron construction in some detail. We 
shall consider also the generalizations of the solutions due to the dilaton, A- 
terms, as well as the solutions for higher gauge groups and in the case of 
axial symmetry. Although very interesting in their own right, these exhibit 
essentially the same features as the BK solutions. 

3.1 Bartnik-McKinnon solutions 

These solutions are known numerically [19,215,55], their existence was estab- 
lished in [293,287,55]. The gauge field potential and the spacetime metric are 
chosen in the form (2.49), (2.50), the field equations are given by (2.51)-(2.54). 
The electric field amplitude a should be set to zero, since otherwise only the 
embedded Abelian solutions are possible [133,121,40]. As a result, the metric 
function a can be eliminated from (2.52) and (2.53) [268] and there remain 
only two independent equations 




2m (w 2 — 1) 2 \ w' w(w 2 — 1) 




(3.2) 



m' = Nw' 2 + 




(3.3) 



2r 2 



Given a solution to these, a is obtained from 



a = exp — 2 / dr — 




(3.4) 
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Fig. 1. w(r) and the effective charge function Q 2 (r) for the lowest BK solutions. 

Eqs.(3.2), (3.3) have singular points at the origin, r = 0, at infinity, r = oo, and 
at r = rh, where N(rh) = 0. Assuming that iV vanishes nowhere (no horizons), 
consider the local power-series solutions to Eqs.(3.2), (3.3) in the vicinity of 
the singular points. Requiring that curvature is bounded one obtains at the 
origin 



At infinity, using the results of [133,121,40], the only nontrivial possibility are 
configurations that asymptotically approach the magnetically neutral solution 



In these expressions b, a and M are free parameters, M being the ADM mass. 
Applying the standard methods one can show that the Taylor expansions in 
(3.5) and (3.6) have a nonzero convergence radius [215,55]. Inserting these 
expressions into (3.4), one obtains 




(3.5) 



(2.55): 




(3.6) 





(3.7) 



at the origin and at infinity, respectively. 
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Fig. 2. Metric functions m(r) and a(r) for the lowest BK solutions. 

The next step is to integrate Eqs.(3.2) and (3.3) with the boundary conditions 
(3.5) and (3.6). The strategy is to extend the asymptotics numerically to the 
intermediate region, where the matching conditions are imposed at some point 
f [266]. There are three such conditions for w, w' and m. This agrees with 
the number of the free parameters in (3.5) and (3.6). It turns out that the 
matching can be fulfilled, provided that the values of the parameters b, a and 
M are restricted to a one-parameter family {b n , a n , M n }, where n is a positive 
integer. This gives an infinite family of globally regular solutions in the interval 
< r < oo labeled by n = 1,2,.... 

The index n has the meaning of the node number for the amplitude w: for 
the n-th solution w(r) has n nodes in the interval r e [0,oo), such that 
w(oo) = (—1)™. For all solutions w is bounded within the strip \w (r)| < 1 (see 
Fig.l). This can be understood as follows: Suppose that w leaves the strip, 
such that w(r) > 1 for some r > (notice that the equations are invariant 
under w — > — w). Then, taking the boundary conditions into account, w must 
develop a maximum outside the strip: w'(r m ) = 0, w(r m ) > 1, w"(r m ) < for 
some r m > 0. However, Eq.(3.2) implies that in the region w > 1 an extremum 
can only be a minimum, which shows that the condition \w\ < 1 cannot be 
violated. 

The behaviour of the metric functions m and o is similar for all n's: they 
increase monotonically with growing r from m(0) =0 and cr(0) = a n at the 
origin to m(oo) = M n and cr(oo) = 1 at infinity, respectively (see Fig.2). Here 
M n is the ADM mass of the solutions and a n is the time deceleration factor 
at the origin. 

Qualitatively the solutions show three distinct regions with the two transition 
zones; see Fig.l - Fig. 3. This distinction is due to the characteristic behavior 
of the effective charge function Q(r) defined by the relation [19] (see Fig.l) 

(3 ,) 
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Fig. 3. On the left: the radial energy density r 2 Tg for the lowest BK solutions. On 
the right: the equations of state T r r /T ° and T^/T °. 

a) In the interior region, r < 1, one has Q(r) 0, the functions w, m and a 
are approximately constant, and the metric becomes flat as r — > 0. The stress 
tensor T^ 4 is isotropic and corresponds to the equation of state p = 3p with 
p oc r 2 ; note that this has the SO(4)-symmetry. The transition zone about 
r = 1 is marked by the high energy density. The metric functions grow rapidly 
in this zone approaching their asymptotic values, while Q(r) reaches the unit 
value. 

b) In the near-field region, 1 < r < R , one has Q(r) m 1, the gauge field is 
well- approximated by that for the Dirac magnetic monopole, w = 0, such that 
p = —p r = p$ with p oc r -4 . The metric is close to the extreme RN metric. In 
the "charge-shielding" transition zone, r ~ R , the effective charge decays to 
zero. 



c) Finally, in the far-field region, r > Rq, the solutions are approximately 
Schwarzschild. The gauge field strength decreases as I/r 3 , the components of 
the stress tensor are p = —p r = 3p$/2 with p oc r~ 6 . The polynomial decay of 
the functions reflects the massless nature of the fields in the problem. 

The above distinction becomes more and more pronounced with growing n. 
The ADM mass M n increases with n from the "ground state" value Mi = 0.828 
to Moo = 1, rapidly converging to the upper limit (see Table 1). The dimension- 
ful mass is 



/ An/Gg 2 M n , and the BK particles are therefore extremely heavy. 
Notice that the discrete energy spectrum emerging in the problem resembles 
that of some quantum system. The analogy becomes even more striking if one 
thinks of the node number n as relating to some wavefunction. The values of 
{b n , a n , M n , a n } for the lowest n are given in Table 1, while those for large n 
can be approximated by 



b n = 0.706 



2.186 e. 



M„ = l 



1.081 e„ 
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Table 1 

Parameters of the BK solutions. 



n 


K 


M n 


a n 


Cn 


1 


0.4537 


0.8286 


0.8933 


0.1264 


2 


0.6517 


0.9713 


8.8638 


0.0208 


3 


0.6970 


0.9953 


58.929 


0.0033 


4 


0.7048 


0.9992 


366.20 


0.0005 


5 


0.7061 


0.9998 


2246.8 


0.00009 




a n = 0.259 e n , 


On = 


0.707 e n 



with e n = exp(— [55]. 

As n — > oo, more and more nodes of w accumulate in the first transition 
zone near r — 1, where N develops a more and more deep minimum, closely 
approaching the zero value. In the near-field zone at the same time, the am- 
plitude and frequency of the oscillations decrease with growing n, while the 
size of the zone stretches exponentially, the whole configurations approaching 
the extreme RN solution: 



W = 0, (7=1, N = 




(3.10) 



However, for finite n the size of the near-field zone is finite, and since w tends 
to ±1 in the far- field region, the solutions are always neutral. 

The limiting solution with n = oo can also be investigated [55,61,294,288]. 
This turns out to be non-asymptotically flat, which can be qualitatively un- 
derstood as follows. The Schwarzschild coordinate system breaks down in this 
limit, and one should use the general parameterization (2.25) with R = R(r) 
(and a = 0). The EYM field equations (see Eqs.(3.24)-(3.26) with A = 0) then 
can be reformulated as a dynamical system. This system admits a separatrix 
that starts from the saddle point with R = 0, w = N = 1, which corresponds 
to the origin, and after infinitely many revolutions ends up at the focal point 
with R — 1, w — N — 0, which corresponds to the degenerate horizon. The 
function R thus changes only in the interval [0, 1]. In contrast, the phase space 
trajectories of solutions with n < oo miss the focal point and propagate fur- 
ther to the the region R > 1, where they approach (3.10) for large n. As 
a result, the configuration that the BK solutions tend to for n — > oo is the 
union of the oscillating solution in the interval R G [0, 1] and the extreme RN 
solution (3.10) for R e [l,oo). Such a configuration cannot be regarded as a 
single solution globally defined for all R > 0. For the oscillating solutions the 
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t =const. hypersurfaces are non-compact and in the limit R — > 1 correspond 
to the cylinders R 1 x S 2 . The four- metric in this limit is the direct product of 
a unit two-sphere and a unit pseudosphere [55]. 

Other solutions to the EYM equations (3.2), (3.3) with the regular boundary 
conditions (3.5) at the origin are also non-asymptotically flat and belong to 
the bag of gold type. This case is in fact generic [55]. The qualitative picture 
is as follows. Using again the general parameterization (2.25) with the gauge 
condition a — 1 the Einstein equation for R(r) is R" = —2w' 2 /R. The regular 
boundary conditions at the origin imply that R(0) = 0, R'(0) = 1. Since 
the second derivative R" is always negative, the first derivative R' is always 
decreasing. As a result, R(r) tends to develop a maximum at some finite r 
where R 1 vanishes. After this the integral curve for R(r) bends down finally 
reaching zero at some finite r, where the curvature diverges. This behaviour is 
generic. Only for the special boundary conditions specified by b = b n in (3.5) 
the integral curves of R(r) bend down slowly enough to be able to escape to 
infinity. For b = b^, R(r) develops an extremum at r = oo, R(oo) = 1. For 
all other values of b in (3.5) R(r) develops a maximum at a finite value of r 
leading to a bag of gold solution. 



3.2 Sphaleron interpretation 

The BK solutions admit an interesting interpretation as EYM sphalerons 
[134,299,143]. This is based on some common features which they share with 
the sphaleron solution of the Weinberg-Salam model [237,208]. In fact, spha- 
lerons exist also in other gauge models [123,43,255] with vacuum periodicity 
[179]. In such theories the potential energy is a periodic function of the wind- 
ing number of the gauge field. The minima of the energy, which are called 
topological vacua, are separated by a potential barrier of a finite height, and 
"sitting" on the top of the barrier there is a classical field configuration called 
sphaleron. Since its energy determines the barrier height, the sphaleron is 
likely to be important for the barrier transition processes when the system in- 
terpolates between distinct vacuum sectors, at least in the electroweak theory 
[217,275]. The winding number of the gauge field changes during such pro- 
cesses, which leads to the fermion number non-conservation due to the axial 
anomaly. 

In what follows we shall introduce the topological vacua in the EYM theory, 
and show that the odd-n BK field configurations can be reduced to either 
of the neighbouring vacua via a continuous sequence of static deformations 
preserving the boundary conditions. The solutions therefore relate to the top 
of the potential barrier between the vacua, which accounts for their sphaleron 
interpretation. 
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3.2.1 Topological vacua 



The sphaleron construction for the BK solutions starts by defining the topo- 
logical vacua in the EYM theory as fields {g^u, A} with zero ADM mass: 
9nv = V^u, A = iXJdXJ -1 . Here r\^ v is Minkowski metric on R 4 and U=U(x l ). 
Imposing the asymptotic condition at spatial infinity [179] 

lim Ufa;*) = 1, (3.11) 



any XJ(x l ) can be viewed as a mappings S* 3 — > SU(2), and the set of all U's 
falls into a countable sequence of disjoint homotopy classes characterized by 
an integer winding number 

k[U] = tr J UdU" 1 A UdU -1 A UdU -1 . (3.12) 

R 3 



As a result, the vacuum fields split into equivalence classes with respect to the 
winding number of the gauge field k[U]. These are called topological vacua. 
Note that the gravitational field is assumed to be topologically trivial. 

A representative of the fc-th vacuum class, {i]^ v , iXJkdU^ 1 } with k[Ujt] = k, 
can be parameterized by 

U fc = exp{ip(r)T r }, where /3(0) = 0, /3(oo) = -2ixk. (3.13) 



The gauge field iU/ c <iU fc 1 is given by Witten's ansatz (2.13) with a = 0, w = 
exp(i/3(r)). 

Note that the asymptotic condition (3.11) imposes the following fall-off re- 
quirement for the components A a of the vacuum gauge field with respect to 
an orthonormal frame: 

A a = o(r _1 ) for r -> oo. (3.14) 



In what follows we shall demand that non-vacuum A's should obey the same 
fall-off condition. 

The vacua with different k cannot be continuously deformed into one another 
within the class of vacuum fields subject to (3.14). However, it is possible to 
join them through a continuous sequence of non- vacuum fields, {g^[A], -A[A]}, 
where the gauge field A[X] obeys (3.14) for all values of the parameter A, while 
<7/«/[A] is an asymptotically flat metric on R 4 . The crucial point is that there are 
such interpolating sequences that pass through the BK field configurations. 
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3.2.2 The interpolating sequence 
Consider a family of static fields [134,319] 



ds 2 = al{r)N x {r)dt 2 - j^dr 2 -r 2 {d$ 2 + sin 2 $d V 2 ), (3.15) 
^[^^U+dU^ + ^UJU: 1 , (3.16) 

where N\(r) = 1 — 2m\(r)/r, and A G [0, ir]. The functions m\(r), cr\(r), and 
U± are given by 



. . sin 2 X f f , 2 . 2 x ( w 1) \ / x , 
m x {r) = — — \w +sm A — ) a x {r)dr, 

°~x( r ) = exp ^—2 sin 2 A J drw^/rj , 

U± = exp{iA(w± l)T r }, (3.17) 

where w coincides with the gauge amplitude for the n-th BK solution. The 
boundary conditions for w ensure that, for any n, 

A a [X] = 0(r) as r -> 0, A a [X] = 0(r" 2 ) as r -> oo. (3.18) 



The metric (3.15) becomes flat for A = 0,7r. The gauge field (3.16) vanishes 
for A = 0, whereas for A = n it can be represented as 

A[n] =iXJdXJ-\ with U = exp{i7r(w - l)T r }. (3.19) 



Now, one has w(0) = 1, w(oo) = (—1)™- Comparing with (3.13) one can see 
that for odd values of n the pure gauge field (3.19) has unit winding number, 
whereas the winding number is zero if n is even. The field sequence (3.15), 
(3.16) therefore interpolates between distinct topological vacua for n odd, and 
between different representatives inside the same vacuum sector for n even. 



3.2.3 The Chern-Simons number 

It is instructive to rederive the above result in a different way [208,134,253,319]. 
Consider an adiabatic time evolution along the family (3. 15), (3. 16) by letting 
the parameter A depend on time in such a way that A(— oo) = 0, A(oo) = ir. 
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Introduce the Chern- Simons current, 



Aivafi 



8tt 2 



-9 



tiA u (v a Ap-jA a A^J , 



whose divergence V CT K a = (l/16n 2 ) tr *F fMU F^ v . For any X(t) in (3.16) the 
temporal component of the gauge field is exactly zero implying that all spatial 
components of vanish. This gives 



t 



(3.20) 



The quantity on the left is iV cs [t] — iV cs [— oo], where iV cs is the Chern-Simons 
number of the gauge field. Since for a pure gauge field N C s coincides with the 
winding number, one has iV"cs[— oo] = 0. Substituting (3.16) into (3.20) and 
evaluating the integral on the right one finally arrives at 



t 

3 



JV cs [A(t)] = — J dt A sin 2 A J dr w'(w 2 - 1) 

-oo 
l-(-l) 



2tt 



(A - sin A cos A). (3.21) 



This shows that for odd values of n the Chern-Simons number changes from 
zero to one as A increases from to n. If n is even then one has iV cs = for 
any A, which corresponds to a sequence of fields within one vacuum sector. 

Now, for A = 7r/2 the fields (3.15)-(3.17) exactly correspond to those of the 
BK solutions. The gauge field (3.15) then can be obtained from the BK field 
via a gauge rotation with U = exp{iir(w — l)T r /2}. The metric function 
m\{r) for A = rr/2 satisfies the Einstein equation (3.3) with the appropriate 
boundary condition and hence coincides with the BK mass amplitude m(r), 
while <J\(r) coincides with a in (3.4). As a result, the configuration (3.15)- 
(3.17) for A = 7r/2 is the gauge rotation of the BK field. 

Summarizing, there are sequences of static fields that interpolate between 
distinct vacuum sectors passing through the BK configurations with odd n. 
Using (3.21), one can determine the Chern-Simons number for these solutions: 
iV cs = 1/2. This shows that the odd-n BK solutions reside in between the two 
vacua and can be reduced to each of them through a continuous sequence 
of static deformations preserving the boundary conditions. This accounts for 
their interpretation as sphalerons. Note also that these solutions admit a zero 
energy fermion bound state [143,319,74]. Moreover, the spectrum of the Dirac 
operator for fermions in the background fields (3.15)-(3.17) for odd values of n 
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exhibits the standard level crossing structure [319]. As the parameter A varies 
from zero to tt, the lowest positive energy level crosses zero at A = n/2 and 
finally replaces the highest energy level of the Dirac see. The fermion number 
therefore changes. For even values of n the BK solutions are topologically 
trivial, since (3.21) gives in this case iV C s = 0. 

3.2.4 The energy and the negative modes 

The above considerations show that the odd-n BK solutions relate to saddle 
points on the potential barrier separating the neighbouring topological vacua. 
For static off-shell fields in GR the potential energy is the ADM mass - pro- 
vided that the configuration is asymptotically flat and fulfills the initial value 
constraints. Thus, for example, for a spherically symmetric, globally regular 
static field configuration the energy is 



Here the second equality on the right is due to the initial value constraint 
Gq = 2Tq (in the dimensionless units chosen), since one has G® = 2m' /r 2 . It 
is worth noting that for any value of A in (3.15)— (3.17) the fields do satisfy 
the initial value constraints, from which the only non-trivial one is Gq = 2Tq. 
In addition, the G r r = 2T£ Einstein equation is also fulfilled. As a result, the 
energy along the field sequence (3.15)-(3.17) is obtained by simply taking the 
limit r — > 00 in m\{r) in (3.17). The result is given by Eq.(3.23) below (with 
(3 — 1). The energy vanishes for A = 0, n and reaches the maximum at the 
sphaleron position A = rr/2 (see Fig. 4). 

The last remark suggests that the BK solutions are unstable. The correspond- 
ing negative mode is given by the derivative of (3.15)— (3.17) with respect to A 
at A = 7r/2. The two facts are essential: a) the energy is maximal at A = ir/2; 
b) the boundary conditions at infinity (3.18) hold for any A, which ensures that 
the negative mode is normalizable. The stability analysis presented in Sec. 5 
below gives the following result: the n-th BK solution has 2n negative eigen- 
modes in the spherically symmetric perturbation sector, of which n modes are 
parity-even and n are parity-odd [323]. The described above negative mode 
belongs to the odd-parity sector. 

3.2.5 The physical picture 

The fact that the BK solutions have more then one negative mode reflects the 
important difference between the EYM and the YMH sphalerons. Sphalerons 
are usually associated with non-contractible loops in the configuration space 



00 




(3.22) 
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Fig. 4. On the left: The ADM mass along the field sequence (3.15)— (3.17) through 
the n = 1 BK solution. On the right: The potential energy (3.23) in the vicinity of 
this solution. 

[237]. Such loops can be obtained from the interpolating sequences considered 
above by identifying the end points. The maximal value of energy along a loop, 
E ma , x , is minimized over all loops to obtain E s = inf {E ma , x }. If a) E s exists and 
is positive; b) there is a loop whose E max is equal to E s , then there is a saddle 
point solution called sphaleron [237]. By construction, this has only one neg- 
ative mode. In reality, however, it is very difficult to show that the conditions 
a) and b) above hold, and this has never been done. At the same time, these 
conditions are essential, since the mere existence of non-contractible loops does 
not imply the existence of any non-trivial solutions [237]. For example, in flat 
spacetime pure YM theory there is vacuum periodicity but there are no static 
finite energy solutions at all. For these reasons a somewhat weaker definition is 
adopted in this text. According to this, sphalerons are only required to relate 
to the top of the potential barrier between the vacua. As a result, they are 
still topologically non-trivial, but can have more that one negative mode, and 
their existence is not guaranteed by the minimax argument. 

The following physical picture seems plausible: The vacuum-to-vacuum transi- 
tions in the EYM theory are suppressed at low energies, whereas the suppres- 
sion is removed at the energy of the order of mass of the n = 1 BK solution. In 
order to justify this picture, one should take into account the fact that, since 
the n = 1 BK solution has more than one negative mode, its energy does not 
determine the minimal height of the potential barrier between the vacua. In 
fact, the latter is zero. In order to see this, it is instructive to consider the 
interpolating sequence (3.15)-(3.17) through the n = 1 BK solution, and to 
generalize it by replacing w(r) by w(flr), where (3 is a parameter. Since the 
boundary conditions do not change, the new sequence still interpolates be- 
tween the distinct vacua, but, unless /3 — 1, no longer passes through the BK 
configuration. For any j3 the fields fulfill the (00) and (rr) Einstein equations, 
which allows one to define the energy as m\(oo), where m\{r) is given by 
(3.17) with w(r) replaced by w(/3r). This gives 
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M((3, A) = /3sin 2 A j L' 2 + sin 2 \^^\ 
o ^ ' 

x exp (^-2(3 2 sin 2 A J ? w ,2 ^j dr, (3.23) 

where w = w(r) corresponds to the n = 1 BK amplitude (see Fig.4). For a 
fixed value of f3 the energy is maximal at A = it/ 2 and vanishes as f3 — > 0, oo, 
so that the minimal barrier height is zero [318,316]. 

In view of the last remark, one can wonder as to why the barrier transitions in 
the EYM theory should be suppressed at low energies [318,316]. Note, however, 
that in the low energy limit the EYM theory reduces to pure YM theory, for 
which transitions between distinct vacua are known to be strongly suppressed, 
despite the fact that inf{E max } = 0.Q On the other hand, the existence of the 
energy scale provided by the sphaleron mass suggests that the transitions are 
unsuppressed above a certain energy threshold. However, the evaluation of the 
path integral is necessary in order to make any definite statements. 



3.3 Solutions with A-term 



A natural generalization of the BK solutions is provided by including the 
A-term into the EYM equations [328,67,307]. It is convenient from the very 
beginning to use the general parameterization of the metric (2.25) with a — 1 
and a = 0. The static EYM-A equations then read 



R"—-^, (3.24) 

(NR)'R r = 2 i^Nw' 2 - (W ' 2 ~ 2 1)2 ) + 1 - A# 2 , (3.25) 

(AV)' (3.26) 
There is also a first integral for these equations 

(NR)'R' = 2 (nw' 2 - ^"/H + 1 - A-R 2 . (3.27) 



5 At finite temperature the total transition amplitude in the YM theory can be 
large due to the high number of thermal quanta participating in the transitions, 
although for each individual quantum the amplitude is small. 
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Fig. 5. On the left: amplitudes w, N, m, and a for the asymptotically de Sitter 
solution with n = 3 and A = 0.0003. On the right, the conformal diagram for this 
solution (the black regions should not be confused with spacetime singularities). 
The definition of the null coordinates U and V is given in [328] . 

The equations admit, for example, the de Sitter solution: w = ±1, R = r, 
N = 1 — Ar 2 /3. For A = the BK configurations fulfill these equations. The 
ra-th BK soliton has the typical size Rq(ti), for R> R the metric being almost 
vacuum. It follows that if the A-term is non-zero but very small, A <C 1/-R 2 ,, 
then the contribution AR 2 to the energy density is negligible. For R < R 
one therefore expects that the solutions do not considerably deviate from the 
BK configurations. In the region R > Rq, however, the effect of A becomes 
significant, which suggests that the metric approaches the de Sitter metric. 
Hence - for sufficiently small values of the cosmological constant - the solutions 
are expected to resemble the regular BK solitons surrounded by a cosmological 
horizon at R ~ 1/ y/X and approach the de Sitter geometry in the asymptotic 



The numerical analysis confirms these expectations [328]. Starting from the 
regular boundary conditions at the origin 



the numerical procedure breaks down at some r h > where N vanishes, which 
corresponds to the cosmological horizon. The local power-series solution at the 
horizon contains four parameters: r h , R h , R' h , and Wh- 



with x = r—r h . Here Wh and N' h are expressed in terms of Rh, R' h , and Wh [328]. 
The strategy is to extend numerically the local solutions (3.28) and (3.29) and 



region. 



w = l-br 2 + 0(r 4 ), R = r + 0(r 5 ), 
N = l + (46 2 -A/3)r 2 + 0(r 4 ), 



(3.28) 



w = w h + w' h x + 0(x 2 ), R = R h + R' h x + 0(x 2 ), 
N = N' h x + 0(x 2 ), 



(3.29) 



32 



impose the matching conditions at some point f, where < f < r^. For any 
A C 1, this leads to a family of solutions in the interval < r < r^, which 
are parameterized by the node number n = 1, 2, . . .. For the n-th solution w 
oscillates n times for < r < r h (see Fig. 5). 

Since the boundary conditions at the cosmological horizon are already fixed, 
the next step is merely to integrate from the horizon outwards to extend the 
solutions to the asymptotic region. For r — > oo one finds 

w = w O0 + ^- + 0(R' 2 ), R = R' oo r + 0{l), (3.30) 

N = 1 ~ 2 -f-^ R2 + ((«4 - 1) 2 - y ° 2 ) h + °( R ^ (3- 31 ) 

where the parameters w^, a, R'^, and M are determined numerically. It turns 
out that Woo 7^ ±1, which means that the magnetic charge does not vanish 
and in the string gauge (2.17) is given by 

P = ^fF=(wl-l)T 3 . (3.32) 

s 2 



The geometry in the asymptotic region is RN-de Sitter. It is interesting that 
the charge determined by the asymptotic behaviour of N does not coincide 
with the one in (3.32). Some numbers are: for the n — 1 solution with A = 0.01 
one has b = 0.452, R h = 16.431, w h = -0.944, Woo = -1.001, M = 0.821 
[328]. The spacetime is topologically non-trivial, and the conformal diagram 
is qualitatively identical to that for the de Sitter solution (see Fig. 5). Similar 
to the de Sitter case, the diagram contains two boundaries R — 0, which now 
correspond to a pair of BK solutions. One can think of the spacetime manifold 
as the de Sitter hyperboloid slightly deformed by masses of two BK particles 
placed at the opposite sides of the spatial section. 

Let us consider now what happens for large values of A. When A grows, the co- 
efficient R'qq in (3.30) decreases until it vanishes for some critical value A crit (n). 
One has A crit (l) = 0.3305. For large r the configurations then approach the 
Nariai solution: w = ±1, R = l/VX, N = -Ar 2 . For A > A crit , R(r) is 
no longer monotone: it reaches a maximum at some finite r e and then starts 
decreasing until it vanishes at some r sing , where N diverges. Such solutions 
belong to the bag of gold type, the 'bag' containing an event horizon and a 
spacetime singularity. If A continues to increase, the position of the maximum 
of R moves towards the origin, while the horizon shifts towards the singularity. 
Finally, for A = A reg (n), the horizon merges with the singularity: — > r sing . 
What remains is a completely regular manifold whose spatial sections are 
topologically S* 3 and have the reflectional symmetry with respect to r e . One 
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Fig. 6. On the left: Change of the topology for the EYM-A solutions. The n = 1 so- 
lution for A = Ai = 0.3304 is asymptotically de Sitter. The one for A = A2 = 0.3305 
is of the Nariai type, while the one for A = A3 = 0.3306 is of the bag of gold type. 
On the right: the regular compact solution for n = 3. 

has A reg (l) = 3/4, in which case the solution is known analytically [109] and 
describes the static Einstein universe 

N=l, R = y/2sm(r/y/2), w = cos^/v^), (3.33) 

which corresponds to the standard metric on R x S 3 . For solutions with n > 1 
one has 3/4 > A rcg (n) > A rcg (oo) = 1/4, and the spatial sections are squashed 
three-spheres. 

No regular solutions exist for A > A rcg (n). In this case N is everywhere positive 
and diverges at r S i ng which is the position of the second zero of R. All EYM-A 
solutions described above are unstable [67]. 

3.4 Stringy generalizations 

It is natural to wonder whether the BK solutions play any role in string theory. 
The most popular stringy generalization of the EYM theory is the EYM- 
dilaton theory (2.30). One can also consider a more general model with the 
(4D) Lagrangian 

£ = ~J + \&<P? - \ ^ (F^Fr - PRgb) - U{4>). (3.34) 

Here R GB = R tlvXr R llvXT - AR^R^ + R 2 is the Gauss-Bonnet term, (3 is a 
parameter, and U((f>) is the dilaton potential. For 7 = 1, U((f>) = 0, (3 = 1 
this corresponds to the toroidal compactification of the low-energy heterotic 
string effective action. The axion can be set to zero provided that the gauge 
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field is purely magnetic. If 7 = 1, (5 — 0, and U(<f>) = (— 1/8) exp(— 20), then 
the theory (3.34) can be obtained via compactification of ten-dimensional 
supergravity on the group manifold [83]. 

In the static, spherically symmetric case the 2D matter field Lagrangian fol- 
lowing from (3.34) is 

L m = L YMD + 2 7 /3e 2 ^ ^ {Na 2 )'(l - NR' 2 ) - R 2 U{<j)), (3.35) 

cr 



where L Y md is given by Eq.(2.31). If f3 — 0, the field equations in the gauge 
R = r are: 



(e^aNwJ = e 2 ^a W ^ \ X \ (3.36) 
{r 2 aN<f>')' = 2 7 ae 2 ^ (nw' 2 + ^{w 2 - l) 2 ) - r 2 aU' (</>), (3.37) 



m' = r l N(f>' 2 + e 2 ^ (nw' 2 + ^(w 2 - l) 2 ) + r 2 f/(0), (3.38) 



(lna)' = r0 ,2 + ^<V 2 . (3.39) 
We shall consider these equations for a number of special cases. 



3.4.1 EYMD theory 

Let us choose U(<f) =0 and (3 = 0. Eqs.(3.36)-(3.39) exhibit then the global 
symmetry 

0^0 + o , r^rexp(70 o ), (3.40) 



which allows one to set 0(0) = 0. The local asymptotic solutions to (3.36)- 
(3.39) are given by 

w = l_6r 2 + 0(r 4 ), m = 0(r 3 ), = 2 7 6 2 r 2 + 0(r 4 ); (3.41) 
w = ± - +0(r~ 2 ), m = M + 0(r- 1 ), 

= 0oo + ^ + O(r- 2 ), (3.42) 

at the origin and at infinity, respectively. Here D is the dilaton charge. The 
symmetry (3.40) implies the existence of the conserved quantity 

Nar 2 {<y In Na 2 - 2(f)}' = C, (3.43) 
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where C is an integration constant. In the regular case one has C — 0, which 
leads to a first integral of the field equations: 



5-00 = exp ( ^ ((f) - ^oo) J , (3.44) 



where goo = ^Vcr 2 . This implies the relation between the dilaton charge and 
the ADM mass [114]: 



D = ~fM. (3.45) 



It follows also that for the string theory value, 7 = 1, the string metric ds 
e 2<t> ds 2 is synchronous. 



Numerical integration of Eqs.(3.36)-(3.39) with the boundary conditions spec- 
ified by (3.41), (3.42) shows that analogues of the BK solutions exist for any 
value of 7 [112,219,35]. These EYMD solitons are also labeled by the node 
number noiw, and the behaviour of the metric and the gauge field is qualita- 
tively the same as in the EYM case. The dilaton is a monotone function. In the 
limit where 7^0 the dilaton decouples and one recovers the BK solutions. 
In the opposite limit, 7 — > 00, the gravitational degrees of freedom decouple 
after the rescaling 

7 — > 00, (j) ~~ 4>lli r 7 r - (3.46) 
Eqs.(3.36)-(3.39) then reduce to those of the flat spacetime YM-dilaton model: 

(e 2 V)' = e 2 ^^^, (3.47) 
(rV)' = 2e 2 ^ ^ /2 + (w2 2 ~ 1)2 ). (3.48) 

Remarkably, these also admit regular solutions with the same nodal structure 
as in the gravitating case [218,36]. The dilaton can therefore play a similar role 
like the Higgs field. For all values of 7 the EYMD solutions have 2n negative 
modes in the spherically symmetric perturbation sector. 

Consider the modification of the EYMD solutions due to the Gauss-Bonnet 
term [113]. It turns out that when f3 in (3.34) is non-zero and small, the qualita- 
tive structure of the solutions does not change. However, if f3 exceeds some crit- 
ical value P CT (n), the solutions cease to exist. One has (5 CT (n) < /3 cr (l) = 0.37. 
As a result, for the string theory value, j3 — 1, there are no regular particle-like 
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solutions. For (3 7^ the scale symmetry (3.40) still exists leading to the con- 
served quantity analogous to that in (3.43). For regular solutions this implies 
the equality M = 7D, but the relation (3.44) no longer holds. 



3.4-2 Gauged supergravity 



For 7 = 1, U((f>) = -(l/8)exp(-20), and (5 = Eqs.(3.36)-(3.39) are inte- 
grable [82,83]. The action (3.34) corresponds then to the truncation of the 
N = A SU(2)xSU(2) gauged supergravity. As a result, one can use the su- 
persymmetry tools to derive the following system of first order Bogomol'nyi 
equations: 



Na 2 = e 2((p -^ ) , (3.49) 
1 + w 2 t M (w 2 - l) 2 i r 2 

rc 



N = ^ + e »^_^ + . (3 . 50) 



'^m^i 1 -^^^ ■ (3 - 51) 



W = -2»^*(l+2e 2 *^i), (3.52) 

which are compatible with the field equations (3.36)-(3.39). The Bogomol'nyi 
equations are completely integrable and admit a family of globally regular 
solutions which can be represented in the form 

dg2 = ^ Sinhp r 2 _ 2 _ R 2 {()){M 2 + ^2 d 2)1 ) (3>53) 

R(p) >■ 

i? 2 (p) = 2pcothp--^-l, w = ±-?-, e 2 * = a 2 ^ (3.54) 
smh p smh p 2 i?(p) 



Here p G [0, 00), and the parameter a reflects the presence of the scale sym- 
metry (3.40). These solutions are stable, preserve 1/4 of the supersymmetries, 
and have unit magnetic charge. When expressed in the Schwarzschild coordi- 
nates, the asymptotics of the solutions are 



N = 1 + + 0(r 4 ), Na 2 = 2e 2 * = a 2 + ^- + 0(r 4 ), 



ro = 1 __ + (A 



iVoclnr, Na 2 = 2e 2(p oc 



4 lnr 



41nr' 



w oc 



(3.55) 
(3.56) 
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at the origin and at infinity, respectively. The geometry is flat at the origin, but 
it is not asymptotically flat. The spacetime manifold is geodesically complete 
and globally hyperbolic. 



3. 5 Higher rank groups 



The simplest generalization of the SU(2) EYM(D) theory to higher gauge 
groups is the SU(2)xU(l) theory [136,111]. However, unless the U(l) field 
vanishes, this does not admit regular particle-like solutions, although black 
holes are possible. The situation improves for larger groups. The SU(N) case 
has been studied most of all. The existence of different embeddings of SU(2) 
into SU(N) leads to several inequivalent expressions for spherically symmetric 
gauge fields. Most of these correspond to models with gauge groups being 
subgroups of SU(N), and there is only one genuinely SU(N) ansatz [333,213]. 
In the static, purely magnetic case this is given by (2.20). The EYM field 
equations read [214] 



r 2 Nw'J + 2(m - rV)w) = - q 0+l )w 3 , (3.57) 

m' = Ng + V, (3.58) 
(\na)' = 2Q/r, (3.59) 

where j = 1, . . . , N-l and % = j(N-j)w]-(J-l)(N-j + l)w]_ 1 +2j-N-l 
with w = wn = 0. One has 

N-l i Af 

G=Y,J(N-J)w?, V = —Y.i- (3.60) 

J=l 3=1 



Solutions to these equations for iV = 3,4 have been studied numerically 
[214,204,206,205,207,194,295]; see [240] for the corresponding existence proof. 
Since the results do not depend considerably on whether the dilaton is in- 
cluded or not, we shall discuss only the pure EYM case. Consider the SU(3) 
theory, when there are two independent YM amplitudes W\ and w 2 . The local 
power-series solutions read 



Wl = \-b ir 2 + b 2 r 3 + 0(r A ), 

W2 = 1 -b x r 2 -6 2 r 3 + 0(r 4 ), m = 0(r 3 ); (3.61) 
Wl = ±(l-^)+0(r- 2 ), 

w 2 = ±(l-—)+0(r- 2 ), m = M + 0(r- i ) (3.62) 
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at the origin and at infinity, respectively. Notice the appearance of the cubic 
terms in the expansions of the w/s at the origin. The numerical matching 
with five shooting parameters, bj, aj and M, gives global solutions. For these 
the behaviour of m and a is similar to that of the BK case, while the gauge 
field amplitudes obey the condition \wj\ < 1 and are characterized by two 
independent node numbers n\ and ni- If w\ = w 2 = w , the equations reduce 
to those of the SU(2) case - up to the rescaling r — > 2r, m — > 2m. As a 
result, there are solutions with the node structure (n, n), corresponding to the 
embedded BK solutions. There are, however, solutions with the same node 
structure but with w 1 7^ w 2 , and they are slightly heavier. Finally, there are 
solutions with arbitrary (ni,n2), where n\ 7^ n 2 - The novel feature is that one 
of the numbers can vanish, while the corresponding amplitude is not constant. 
The lowest value of the mass M(rii,n 2 ) is M(1,0) = 2 x 0.653, then comes 
M(2,0) = 2 x 0.811, followed by the doubled mass of the embedded n = 1 
BK solution, M(l, 1) = 2 x 0.828. The next is the mass of the genuine SU(3) 
solution, M(l, 1) = 2 x 0.847, and so on. In the limit n — > 00 the (0,n) and 
(n,n) sequences display the oscillating behavior in the inner region, outside of 
which the geometry tends to the extreme RN solution with magnetic charges 
P = \/3 and P = 2, respectively. For comparison, in the SU(2) case the 
limiting solution is characterized by P = 1. If the dilaton is included, then the 
limiting solution in the exterior region corresponds to the extreme dilatonic 
black holes with magnetic charge P. 

In the SU(4) case the nodal structure is determined by (711,712,713) [194]. The 
diagonal sequence (n,n,n) contains three different types of solutions. The first 
one, with w± = w 2 = W3, corresponds to the n-th BK solution rescaled by the 
factor vTO. For solutions of the second type two of the three w/s are equal, 
while for those of the third type all the w/s are different. For n — > 00 these 
solutions approach in the exterior region the extreme black hole configuration 
with the charge P = y/lO. Other sequences tend in a similar way to the 
extreme solutions with other charges, whose values can be classified using the 
algebraic structure of SU(4). A similar classification can also be carried out 
in the general SU(N) case [207]. One should have in mind, however, that for 
finite n the regular solutions are always neutral. For n — > 00 they tend to 
a union of two solutions, one of which is non-asymptotically flat and exists 
in the interval < r < P, while another one is the extreme RN solution in 
the region r > P. For the regular SU(N) solution with dilaton the equality 
D = 7M still holds. 

All known regular EYM solutions for higher gauge groups are unstable [71,73]. 
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3.6 Axially symmetric solutions 

The EYM(D) field equations in the axially symmetric case are obtained by 
using the ansatz (2.22) for the gauge field and that in (2.21) for the metric. 
These are rather complicated partial differential equations for the seven real 
amplitudes in (2.21), (2.22). Due to the residual gauge invariance of (2.22), the 
number of independent amplitudes is six. The equations admit interesting so- 
lutions generalizing the BK particles to higher winding numbers [199,195,201]. 
Similar generalizations exist also in flat space for monopoles [269] and spha- 
lerons [196,197]. The basic idea can be illustrated as follows. The axial ansatz 
(2.22) for the gauge field covers, in particular, the spherically symmetric case. 
Specifically, choosing in (2.22) v — 1 and 



with z + ip = re 1 ^ and w(r) = w*(r) and omitting the tilde sign, the expres- 
sion reduces to the spherically symmetric ansatz (2.13) with a = a r = 0. In 
a similar way, the axially symmetric line element (2.21) can be used in the 
spherically symmetric case. As a result, given a static, spherically symmetric 
solution one obtains a solution for the axial EYM equations with v — 1. The 
next step is to use this solution as the starting point in the numerical iteration 
scheme with v — 1 + 5u. The iterative field configurations must satisfy cer- 
tain regularity conditions at the symmetry axis and at infinity. The iterations 
converge, and repeating the procedure one obtains in this way solutions for 
arbitrary v. The physical values of v are integer. This finally gives the gener- 
alized BK solutions characterized by a pair of integers (n, u), where n is the 
node number of the amplitude w in (2.22). 

For v 7^ 1 the solutions are not spherically symmetric. The contours of 
equal energy density Tq are 2-torii and squashed 2-spheres. The mass of 
the solutions, M(n,u), increases with v. For example, for n — 1 one has 
M(l,l) = 0.828 (the BK solution), M(l,2) = 1.385, M(l,3) = 1.870; while 
for n = 2 one finds M(2, 1) = 0.971, M(2, 2) = 1.796, M(2, 3) = 2.527 and so 
on [201]. If the winding number v is fixed while the node number n tends to 
infinity, then the mass approaches the value M(oo, v) = v. In this limit the 
solutions exhibit a complicated oscillating behaviour in the interior region, 
outside of which region the configurations approach the extreme RN solution 
with the magnetic charge v. 

One can think of the solutions with higher winding numbers as describing non- 
linear superpositions of v BK particles aligned along the symmetry axis. This 




a p (p,z) + ia z (p,z) 



w(r) — 1 



(3.63) 



r 
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is supported by the fact that their Chern-Simons number, which is computed 
in the same way as for the BK solitons, is equal to u/2. The solutions can be 
generalized by including a dilaton field for an arbitrary value of the coupling 
constant 7. It turns out that the metric-dilaton relation (3.44) holds in the 
axially symmetric case too. This implies the relation (3.45) between mass and 
the dilaton charge. For 7 — > 00 gravity switches off, but the solutions survive 
in this limit [198]. Although the corresponding stability analysis has not been 
carried out yet, it is very likely that all these EYMD solitons are unstable. 

It remains unclear whether the procedure described above gives all solutions 
in the static, purely magnetic, axially symmetric case. It is possible that con- 
figurations violating the circularity condition can exist as well. In particular, 
it is unclear whether the BK solutions exhaust all possibilities for v — 1. 



4 Non-Abelian black holes 

The intuitive physical idea behind the original no-hair conjecture [277] ap- 
parently was that only exact physical symmetries like gauge symmetries can 
survive in a catastrophic event like a gravitational collapse. Associated with 
gauge fields there are conserved charges, which can be measured using the 
Gauss flux theorem. On the other hand, all physical quantities which are not 
coupled to the corresponding gauge fields, like baryon number, are not strictly 
conserved. As a result, they either disappear during the collapse or become 
unmeasurable. The proof of this conjecture was given for a number of spe- 
cial cases. In such a proof it is not the process of collapse that is usually 
considered but the result of it - a stationary black hole spacetime with cer- 
tain matter fields. It turns out that black holes cannot support linear hair 
for scalar [86,25,303], spinor [157,304] and massive vector [24,23] fields; see 
[149,22,21,242,164,239] for a more recent discussion. A similar proof was given 
also for some non-linear matter models [4,22,242,162,163,298]. 

An interesting confirmation of the no-hair conjecture was found in the fer- 
mionic sector of the N=2 supergravity, where there are black holes with an 
external spin-3/2 gravitino field [6]. The latter, being subject to the Gauss 
law, gives rise to the fermionic supercharge of the black hole, and this does 
not contradict the conditions of the conjecture. Note that in the bosonic sector 
supergravity black holes generically support external scalar fields associated 
with dilaton and moduli fields [138], and these are not subject to a Gauss 
law. A similar example is provided by black holes with conformal scalar field 
[42,26]. However, in all these cases the scalars are completely parameterized 
by black hole mass and gauge charges and hence cannot be considered as inde- 
pendent - sometimes they are called secondary hair. In addition, the Abelian 
supergravity black holes are subject to the uniqueness theorems [62,60]. All 
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this solidified the belief in the universal validity of the no-hair conjecture. 

The first example of the manifest violation of the conjecture was found in the 
EYM theoryQ. In the early study of the EYM system it was observed that 
for any solution (g^, A^) of the Einstein-Maxwell theory there is a solution 
(gnu, °f the EYM theory with the same metric and = T^4 M (x), where T 
belongs to the Lie algebra of the gauge group [340]. Therefore, either directly 
solving the EYM equations or starting from the Kerr-Newman configurations, 
a number of embedded U(l) black hole solutions for the EYM field equations 
had been found [13,87,330,261,184,189-191]. These were called colored black 
holes [261]Q. In the spirit of the no-hair and uniqueness conjectures, it was 
assumed for some time that these solutions exhaust all stationary EYM black 
holes. A partial confirmation of this assumption was given by the "non-Abelian 
baldness theorem" proven for the SU(2) gauge group [133,121,40]: all static 
EYM black hole solutions with finite colour charges are the embedded Abelian 
ones. This assertion, however, left open a possibility to have essentially non- 
Abelian EYM black holes in the neutral sector, and such solutions were found 
in [324,215,33]. For these black holes the only parameter subject to the Gauss 
law is their mass, and for a given value of the mass there can be several 
different solutions. The no-hair conjecture is therefore violated. 

The existence of one example triggered a broad search for other hairy black 
holes in various non-Abelian models. This has led to their discovery inside 
magnetic monopoles [258,222,54,56], Skyrmions [117,38], as well as in a num- 
ber of other systems. For the reader's convenience we list here the families of 
hairy black holes known up to date: 

(1) Black holes with pure Yang-Mills fields. Apart from static and spher- 
ically symmetric SU(2) EYM black holes [324,326,215,33], this family 
includes also static and axially-symmetric solutions [200,202], stationary 
and non-static configurations [327,69], solutions for higher gauge groups 
[136,204,207], and the generalizations including a dilaton and the higher 
curvature terms [112,219,306,111,257,309,302,206,205,194,295,185,187,7]. 

(2) Black holes with Yang-Mills-Higgs fields. These are magnetically charged 
solutions with the triplet Higgs field [258,222,54,56,5], their analogues for 
a more general matter model as well as the generalizations to the non- 
spherically symmetric case [273,271,332], and the neutral solutions with 
the doublet Higgs field [149,231,308,301]. 

(3) Skyrme black holes [228,117,172,38,203]. 



6 Hairy black holes in the Einstein-Skyrme model had been described in [228] be- 
fore the EYM solutions were found, which paper, however, has remained almost 
unknown. 

7 Unfortunately, despite the fact that they are neutral, the same name stuck later 
to the non-Abelian EYM black holes. 
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In view of this variety of hairy black holes one can wonder as to how to reconcile 
their existence with the original arguments behind the no-hair conjecture. 
One possibility, which can probably be applied for unstable solutions, is to 
argue that hairy black holes cannot appear via gravitational collapse. Another 
option, which is suggested by the examples of stable black holes inside solitons, 
is to assume that the topological charge cannot disappear during the collapse. 
This is supported by the fact that in all known examples the size of black holes 
inside solitons is bounded from above, otherwise no hairy solutions exist. It 
seems then that a black hole cannot swallow up a topological object with the 
typical size exceeding that of the black hole. 



4-1 EYM black holes. The exterior region 



The EYM black hole solutions satisfy the same equations as the regular BK 
solutions [324,326,317,215,33]: 

AV< + ( *2 - W - = W{w2 : l \ (4.1) 

I ny ry A J ry ry £1 

m' = Nw» + (4.2) 
The boundary conditions at infinity are still given by (3.6) 

w = ±(l-^+0 (r~ 2 ) , m = M + (r~ 3 ) , (4.3) 

which are now supplemented with the requirement that there is a regular event 
horizon at r = r\ > 0: 

N(r h ) = 0, N(r) > for r > r h . (4.4) 



The regularity assumption implies that all curvature invariants at r = 
are finite. The local power-series solution to (4.1), (4.2) in the vicinity of the 
horizon reads 



w = w h + Wh{W l 1} (r - r h ) + 0((r - r h ) 2 ), (4.5) 

h 

N= ^- h ( 1 - {Wl ~\ ^ " + 0((r " r " )2) ' (4 ' 6) 

It follows from the second condition in (4.4) that r\ > (w\ — l) 2 . 
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Fig. 7. On the left: the amplitudes w, m and a for the non-Abelian black hole 
solutions with = 1. On the right: the effective charge Q 2 and the radial energy 
density t 2 Tq. 



For a given r h > Eqs.(4.1)-(4.6) define a non-linear boundary value problem 
in the interval < r < oo. The following simple argument is in favour of the 
existence of a non-trivial solution at least for ^> 1. Passing everywhere to 
the radial coordinate x = (r — r^/r^, some terms in the equations acquire 
the factor and can be omitted for large rh, if only they are bounded 
for x > 0. As a result, the Einstein equation (4.2) decouples and admits the 
solution m = Th/2 corresponding to the Schwarzschild metric. The remaining 
YM equation (4.1) reads 

( XW ' V = W ( W2 - 1) r4? N 



where the differentiation is with respect to x. The following solution to this 
equation is known [48]: 

2x-l-V3 , in . 

w = ± —=. 4.8 

2^ + 5 + 3^ 



This describes a regular YM hair on the Schwarzschild background. This 
shows that the non-trivial solutions to the full system of equations are likely 
to exist, at least for r h 3> 1. 

In order to obtain the solutions to the full problem (4.1), (4.2) the procedure is 
to numerically extend the asymptotics (4.3), (4.5) to the intermediate region. 
The matching conditions for w, w' and m determine then the three parameters 
M, a and Wh in (4.3), (4.5). As a result, for any given > 0, one finds 
a sequence of the global solutions in the interval r^ < r < oo. These are 
parameterized by the node number n of w. The non-Abelian black holes can 
therefore be labeled by a pair (rh, n), where > 0, n — 1, 2, The existence 
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of these solutions has been established in [292,55]. 



For any (r^, n) the behaviour of the amplitudes w, m and a in the exte- 
rior region is qualitatively similar to that for the regular BK solutions (see 
Fig. 7). The amplitude w starts from some value < w n (r h ) < 1 at the 
horizon and after n oscillations around zero tends asymptotically to (— l) n . 
Considerations similar to those used in the regular case show that one has 
|u>(r)| < 1 everywhere outside the horizon. The metric functions m and a 
increase monotonically with growing r from m(r/j) = r^/2 and cr(r/j) = cr n (rh) 
to m(oo) = M n {r- h ) and cr(oo) = 1, respectively. In the asymptotic region, 
r — > oo, the geometry is Schwarzschild with the mass M = M n (r/ l ) depending 
on r/j and n. Since the YM field strength decays asymptotically as 1/r 3 , the 
Gauss flux integral vanishes and the node parameter n cannot be associated 
with any kind of YM charge of the black hole. 

One can qualitatively distinguish between two regions in the parameter space 
of the solutions: a) r\ > 1, and b) < r\ < 1. This is due to the comparison 
of the "bare" mass of the black hole and its total ADM mass. The bare mass 
is m(r h ) = r h /2. Since m(R) can be thought of as the total energy, including 
the gravitational binding energy, confined in the region r < R, the bare mass 
is the total energy trapped inside the event horizon. Using (4.2), the ADM 
mass M = ra(oo) can be represented as the sum of the bare mass and the 
"dressing" mass 



the second term on the right being the contribution of the matter distributed 
outside the horizon. 

a) Th > 1. In this case the matter term in Eq.(4.9) is small, M n (rh) ~ rh/2, 
and the YM field almost does not influence the geometry outside the horizon. 
The YM field energy increases slightly with growing n, the ADM mass being 
(see Fig. 8) 



Here coincides with the mass of the RN black hole with unit charge. The 
following approximation is very good for n > 2 and r > 1: M n (r h ) ps 
(see Fig. 8). This is due to the fact that for large values of n the amplitude w 
oscillates in a small vicinity of zero. The charge function Q defined by Eq.(3.8) 
is then close to unity, while metric is approximately RN. For n — > oo the metric 
converges pointwise to the RN metric. 



oo 




(4.9) 



M 1 (r h ) < M n (r h ) <M OQ = ^ + ^-. 

2 2r h 



(4.10) 
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Fig. 8. The ADM mass M n (r^) and the temperature screening coefficient P n (fh) 
versus the event horizon radius. For n = solutions are Schwarzschild. 

For r h 3> 1 the matter term in Eq.(4.9) is of the order of l/r h and the back 
reaction of the YM field on the spacetime geometry becomes negligible. The 
solutions in this limit reduce to the bound states of the YM field on the fixed 
Schwarzschild geometry. These are described by Eq.(4.7), whose solution for 
n — 1 is given by Eq.(4.8), while those for n > 1 are known only numerically. 

b) < rh < 1. As rh decreases, the effect of the YM field on the geometry 
becomes more and more pronounced. For ~ 1 the "bare" and "dressing" 
masses are of the same order of magnitude, such that M n (r h ) ~ r h . For so- 
lutions with r h < 1 the role of the YM field is dominant. The solutions then 
look like small black holes dressed in the YM "coat" . As r h approaches zero 
the external field configuration becomes more and more close to that for the 
regular BK soliton with the same value of n. The ADM mass for < 1 varies 
within the following range: 

Mi = 0.828 < M n (r h ) < M 00 (r^) = 1. (4.11) 

Here Mi is the mass of the ground state BK solution (see Fig. 8). For any 
rh < 1, M n (rh) converges to the unit value for n — > oo. In this limit an infinite 
number of nodes of w accumulates near r = 1. As a result, similar to the 
situation in the BK case, the solutions tend for n — > oo to a union of an 
oscillating solution in the interval r h < r < 1 and the extreme RN solution 
for r > 1. 

In the limit — > the event horizon shrinks to zero and the black hole 
solutions converge pointwise for r > to the regular BK configurations. 

The first law of black hole physics for the EYM black holes was considered 
in [170,169], the thermodynamics was discussed in [252,306]. The Hawking 
temperature of the EYM black holes can be computed by analytically contin- 
uing the metric to the imaginary time and requiring the absence of the conical 
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singularity. The result is T = K, sg /2ir, where the surface gravity K sg = <JhN' h /2 
can vanish only for the extreme RN solution [289,290]. Using (4.6) one obtains 

r =wi (412) 

where the second factor on the right is the temperature of the Schwarzschild 
black hole with radius r^, while /3„(r^) = o"(r^)(l — (w\ — l) 2 /r^) (see Fig. 8). 
The presence of I3 n {rh) in this formula leads to the existence of a (narrow) 
region in the parameter space, ri(n) <r h < r 2 (n), for which the temperature 
increases with growing M and hence the specific heat is positive [306]. 

All EYM black holes are unstable and have 2n negative modes in the spher- 
ically symmetric perturbation sector (see Sec. 5 below). Note that for black 
holes the analogy with sphalerons does not exist. The reason is very deep - 
black holes are associated with thermal states and not with excitations over 
pure vacuum states. For black holes one cannot define classical YM vacua even 
formally, since the topology is different and pure gauge fields are no longer 
characterized by integer winding numbers. 

4-2 EYM black holes. The interior structure 

From the conceptual point of view the problem of finding the solutions in the 
interior region r < r\ is simple. One should merely integrate the field equa- 
tions (4.1), (4.2) inwards starting from the event horizon, since the boundary 
conditions (4.5) at are already specified by the exterior solutions. The re- 
sult, however, turns out to be quite bizarre [115,59]. Therefore, before passing 
to the numerical analysis, some preliminary analytical considerations can be 
useful. 



4-2.1 Special solutions 

It turns out that the interior solutions generically do not have an inner horizon. 
Indeed, supposing that such a horizon exists at some r_ < r^, the solution in 
the vicinity of r_ is given by (4.5) with r^ and Wh replaced by r_ and some 
■u;_, respectively. One can extend this solution outwards to match the one 
that propagates inwards from r h . Now, in order to fulfill the three matching 
conditions for w, w' and m one needs at least three free parameters. At the 
same time, only r_ and W- can be used for this, since the solution coming 
from r/j is completely specified by the boundary conditions at rh- One can use 
rh as the third matching parameter. As a result, it might happen that for some 
special values of the event horizon an inner horizon exists too. However, for an 
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arbitrary the matching is impossible. This shows that the non-Abeian black 
holes generically do not possess inner horizons [128,130,58,57] and therefore 
have a spacelike singularity, as suggested by the strong cosmic censorship 
hypothesis. 

Next, one can study the power series solutions to (4.1) and (4.2) in the vicinity 
of the singularity r = 0. One finds three distinct types of local solutions, but 
none of them are generic. First, there are local solutions with a Schwarzschild 
type singularity: 



w = ± ( 1 - br 2 + b 2 A + 0(r 6 ), (4.13) 

V 30m J 

m = m (1 - 4b 2 r 2 + 86 4 r 4 ) + 2b 2 r 3 + 0(r 6 ), (4.14) 

with arbitrary m and b. Since the inner horizon is generically absent, a solu- 
tion of this type should match the one propagating inwards from r^. However, 
the number of free parameters is not enough for this. 

Next, there are solutions with a RN-type singularity: 



w = w + W ° r 2 + cr 3 +0(r 4 ), (4.15) 
2(1 — Wq) 

(w 2 - 1 "l 2 

m = _^o L + mo + (r 2 ). (4.16) 

2r 

Here the number of free parameters is three: m , c, and w ^ ±1, which agrees 
with the number of the matching condition. However, one has iV — > +oo 
for r — * 0, such that iV is positive near the singularity. This requires the 
presence of an inner horizon, which generically does not exist, hence this type 
of singularity is also generically impossible. 

Finally, there is a branch of local solutions containing only one free parameter, 
w [115,59]: 



w = w ± r - - 7° u r 2 + (r 3 ) , (4.17) 
(w 2 - I ) 2 

m= l^o iL ±2wo K 2 -l) + 0(r). (4.18) 

2?" 

For these the singularity is spacelike, and iV oc — (w 2 , — l) 2 /r 2 as r — > 0, which 
formally corresponds to a RN metric with "imaginary charge" . The geometry 
is conformal to R 2 x S* 2 . It is worth noting that solutions of this type were 
found [259] in the study of the phenomenon of mass inflation [262] . 



48 



Summarizing, none of local power-series solutions contain enough free param- 
eters to fulfill the matching conditions. Note also that the possibility of having 
a singularity at a finite value of r can be ruled out [291]. 

To recapitulate, solutions in the interior region do not generically possess in- 
ner horizons and do not exhibit a power law behaviour in the vicinity of the 
singularity. One can wonder then how these solutions look like. The answer 
is provided by a numerical analysis. First of all, it turns out that the special 
solutions described above can indeed be obtained by fine tuning of the param- 
eters - apart from those in (4.17), (4.18). For n — 1 and r h (l) = 0.6138 the 
singularity is Schwarzschild-type, such that there is no inner horizon. For all 
n > 2 there are special solutions which have RN-type singularity with an in- 
ner horizon at some r_(ra). For example for n = 2 one has r^(2) = 1.2737 and 
r_(2) = 0.0217; while for n = 3 one finds r h (3) = 1.0318 and r_(3) = 0.0894. 
It seems that there are also other special solutions with the Schwarzschild-type 
singularity, but these are very difficult to obtain numerically [59]. 



4-2.2 Generic case 

Coming to the generic behaviour of the solutions in the interior region, this 
turns out to be quite unusual. Specifically, choosing an arbitrary external 
solution, with n — — 1, say, and integrating inwards from the event horizon, 
the metric functions m and a and the derivative w' exhibit in the interior region 
violent oscillations, whose amplitude and frequency grow without bounds as 
the system approaches the singularity. At the same time, since everything 
happens at the very short scale, the amplitude w is almost constant in the 
internal region. 

During a typical oscillation cycle the qualitative behaviour of the solution is 
as follows (see Fig. 9). The k-ih cycle (k = 1,2...) starts at the k-ih local 
minimum of the mass function m at some where m(r^) is very close to zero. 
Then, as r decreases from to some Rk < rk, m grows exponentially fast until 
it reaches a very large value at Rk- At the same time a decreases by many 
orders of magnitude. In the region r < R k the mass function reaches a hori- 
zontal plateau where it is almost constant, m ~ M k = m(R k ); similarly for a: 
a PS a k = cr(Rk)- The plateau stretches as far as many orders of magnitude of 
r towards the singularity, during which period the geometry is approximately 
Schwarzschild with a very large mass Mk- Then, at the end of the plateau, m 
catastrophically falls down reaching a very deep minimum at some r fc+1 , after 
which the next oscillation cycle starts. 

The oscillation amplitude for each subsequent cycle is exponentially large com- 
pared to that for the preceding one. As a result, it is extremely difficult to fol- 
low numerically more than the first two or three cycles. It turns out, however, 
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Fig. 9. On the left: the mass function m for the interior solutions with r/j = 1 and 
n = 1,2, 3. The amplitude w is shown for the solution with n = 3. On the right: the 
second oscillation cycle for the n = 3 solution. 

that a simple analytic approximation allows one to qualitatively understand 
the generic behavior of the system for r —>■ [115,59]. Such an approximation 
is based on the truncation of the field equations (4.1) and (4.2) due to the 
numerical observations that for r<r/, one has 



w ~ const. 7^ ±1, 
(w 2 - l) 2 



> 1, 



(w 2 — l) 2 w' . 9 , J . 

i > w{w 2 - 1). (4.19) 



Neglecting in Eqs.(4.1) and (4.2) the small terms compared to the large ones 
the equations reduce to the following dynamical system: 



x = e y - 1, 

y = l+e y -2e 2x , 



whereQ 



e x = w' 



(4.20) 



{w 2 - 1) 
Nr 2 



(4.21) 



and the differentiation is with respect to r = — ln(r/ro), with tq being a 
constant. We are interested in the behaviour of the solutions for r — > oo. 
First of all, let us notice that there is only one critical point of the system, 
(x,y) = (0,0), which is an attractive center for r — ► — oo. As a result, the 

8 The sign of w 1 is definite for r <C r^, and so one can choose w' > utilizing the 
symmetry w — > — to. 
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Fig. 10. The solution to the dynamical system (4.20). 

behaviour of the solutions for large and negative values of r is known: the 
trajectories starting close to the critical point spiral outwards when r increases 
(see Fig. 10). Now, the crucial fact is that this spiraling motion can never 
stop. One can show that, starting from any point on the plain, the trajectory 
cannot escape to infinity for finite r and performs the full revolution around 
the center within finite r. The spiraling behaviour is therefore generic and 
persists forever. Note that the critical point around which the spiraling occurs 
corresponds to w = Wq, N = —(wq — l) 2 /r 2 . For small r this agrees with 
the local solution (4.17), (4.18). Dynamically, however, this solution cannot 
be reached. Indeed, the trajectories approach the critical point only for r = 
— ln(r/r ) — > — oo, such that r — ► +oo, contradicting the assumption that r 
is small and that r < r^. 

The solution to the dynamical system (4.20) shown in Fig. 10 illustrates nicely 
the various phases of a typical oscillation cycle described above. After each 
revolution, the picture in Fig. 10 repeats itself at an exponentially blown up 
scale [1150. The revolutions go on forever and the system approaches the sin- 
gularity via an infinite sequence of more and more violent oscillation cycles. 
At the beginning of each cycle iV is very close to zero, such that the system 
exhibits an infinite sequence of "almost" inner horizons, but the true Cauchy 
horizon never appears. Note that the oscillatory nature of the solution re- 
sembles somewhat the well-known situation in some cosmological models [27]. 
However, the solutions are not chaotic. This follows from the existence of the 
effective description in terms of the two-dimensional dynamical system, since 
chaos cannot occur in two dimensions. Inside black holes, a similar oscillatory 
behaviour was observed in [259] for a homogeneous mass inflation model with 
radial null radiation. 

The values of amplitudes and periods of the oscillations can be estimated 



This opens a possibility for numerically integrating the equations by making a 
logarithmic substitution of the variables after each cycle [344] 
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analytically [131]. The essential quantities for this are Xk = ij-k/Rk) 2 3> 1, 
where and Rk correspond to the beginning of the fc-th cycle and the fc-th 
plateau, respectively. One can show that the x^'s obey the following recurrence 
relation 

x k+1 = -3-, (4.22) 



thus constituting an exponentially diverging sequence. In terms of x k one has 

^ = x k e-*"' 2 , (4.23) 



which can also be viewed as the ratio of the neighboring oscillation periods, 
since r k ^> The increase in M k and decrease in cr fe can be expressed as 



M k e x ^ 2 a k+1 __ Xk/2 



M k -i X k ' cr k 



= e~ Xk/2 , (4.24) 



respectively. One can see that the ratios of the oscillation amplitudes and 
frequencies for the neighbouring cycles are exponentially large, where the ar- 
gument of the exponent grows exponentially for each subsequent cycle. 

Note finally that the exterior structure of the EYM solutions does not change 
considerably after adding some additional matter, like a dilaton or a Higgs 
fields, say. However, the solutions in the interior region change completely 
[131,59,344,127,130,129,280] and exhibit a regular power-law behaviour near 
the singularity. In this sense the oscillatory character of the interior solutions 
described above distinguishes the EYM black holes among all other known 
hairy black holes. Note also that, although the EYM black hole solutions in 
the exterior region are unstable, their interior oscillatory structure is stable 
with respect to non-linear spherically symmetric perturbations [116]. 



4-3 Non-Abelian dilaton black holes 



The EYM non-Abelian black holes can be generalized to include a dilaton 
[112,219,306]. The field equations are then given by (3.36)-(3.39). The bound- 
ary conditions at infinity are specified by (3.42), while at the horizon one 
has 



w = w h + w' h x + 0(x 2 ), (j) = (j) h + <p' h x + 0(x 2 ), N = N' h x + 0(x 2 ) (4.25) 
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with x = r — Thi and w' h , <f>' h , N' h being determined by r^, Wh and fa. In the 
black hole case the scale transformations (3.40) change r^. However, since the 
solutions exist for any r h and fa, one can set fa = without loss of generality. 
The numerical integration gives for any 7 < 00 a family of black hole solutions 
parameterized by (r^,n) as in the 7 = case. The behaviour of w, m, and 
a is qualitatively the same as for the EYM solutions, while the dilaton is a 
monotone function. 

The conservation law following from the scaling symmetry exists in the black 
hole case too, but now the integration constant in (3.43) does not vanish and 
is given by C = r ya} l N' h r'f h . As a result, the relation (3.45) no longer holds, but 
instead one finds 

D = 1 {m-^ K sg A} , (4.26) 

where n sg = ahN' h /2 is the surface gravity and A = knr\ is the area of the 
event horizon. 

The thermodynamics of dilaton black holes depends on value of the dilaton 
coupling constant 7 [306]. For small values of 7 there is a region in the pa- 
rameter space for which the specific heat of the solutions is positive. For large 
7 this region shrinks to zero. The stability behaviour of the solutions, on the 
other hand, do not depend on 7: the number of negative modes is the same 
for all values of 7. 

What is entirely different from the EYM case is the behaviour of solutions in 
the interior region. First of all, for 7 7^ inner horizons cannot exist even after 
fine tuning of the parameters [280]. Indeed, the expression in (3.43) reduces 
at the event horizon to C = ^UhN^r^, which is positive. At the inner 
horizon r_ one has C = 70"_iV"ir 2 , which is negative. This contradicts the 
fact that C is constant, and hence inner horizon cannot exist. The numerical 
analysis reveals the following characteristic picture in the interior region, which 
is generic for all 7 7^ 0: the mass function m, after some oscillations near the 
horizon, exhibits the power-law behaviour all the way down to the singularity 
[127,280]. The power-series solution in the vicinity of the singularity for 7 7^ 
is generic: 

w = w + ar 2(1 " A) , m = fir'^ , = c + In (V A ) , a = o x r x2 , (4.27) 

since it contains the maximal number of free parameters: w , a, c, /i, A, and U\. 
As 7 decreases, the solutions develop more and more oscillation cycles in the 
interior region. However, as long as 7 7^ 0, the oscillations are always replaced 
by the power-law behaviour for r — > 0. This can be illustrated as follows 
[280]. For 7 = the interior oscillating solution can effectively be described 
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by the two-dimensional dynamical system (4.20), whose trajectories evolve in 
the plane (x,y). For 7^0 there is an additional degree of freedom due to 
the dilaton, and the reduced system becomes three-dimensional. As a result, 
after a number of revolutions in the vicinity of the plane, which corresponds 
to several oscillation cycles, the trajectory drifts away towards a critical point 
of the equations that lies outside the plane. This changes the character of the 
solution. 

The EYMD black holes can be generalized in the context of the theory with the 
Gauss-Bonnet term [187,309]. The string theory value of the coupling constant 
in (3.34), /3 — 1, is allowed for black holes. The solutions are parameterized by 
three numbers: r^, n). This is because the event horizon for /3 7^ exists 
only if Th > ?"o > 0, which implies that one cannot use the scale symmetry to 
set <fih = 0. One has r ~ (3 exp(2^(J) h ). For n — 0, when the YM field strength 
vanishes, the solutions describe black holes with dilaton hair [185,186,7]. 



4-4 Higher gauge groups and winding numbers 

The SU(2) non-Abelian black holes admit generalizations to higher gauge 
groups [136,111,114,214,204,206,205,207,195,194,295]. These exist for all val- 
ues of the dilaton coupling constant 7. In the maximal SU(N) case, when all 
the gauge field amplitudes Wj are non-trivial, the solutions are neutral and 
can be parameterized by (rh,rij), where rij is the node number of Wj. The 
regular limit is recovered for r h — > 0. The solution space exhibits bifurcations 
in this limit, since for — > there exist different solutions with the same 
node structure which merge as grows. If some of the n/s tend to infinity, 
the solutions approach the RN (for 7 = 0) or charged dilaton (for 7 7^ 0) black 
hole solutions. 

One of the new features is that when some of the w/s vanish identically, 
the U(l) charges arise. The non-Abelian baldness theorem hence does not 
generalize to higher gauge groups. For example, setting in the SU(3) equations 
(3.57)-(3.59) w 1 = 0, the solutions for w 2 show the usual nodal structure. 
However, since w 1 does not obey the boundary condition in (3.62), there is a 
magnetic charge P = for all n 2 . This is reflecting in the behaviour of N for 
r — > 00: N = 1 — 2M/r + P 2 /r 2 + 0(l/r 3 ). For higher gauge groups one can 
obtain in this way a variety of charged solutions, in particular, solutions with 
an electric charge (for a more general ansatz for the gauge field). For charged 
EYM black holes the event horizon radius is bounded from below, > P, 
since otherwise the singularity is naked. The limit where — > for the 
charged solutions is therefore impossible. Solutions with = P are extreme 
black holes with degenerate horizon. For charged solutions with dilaton the 
r h — > limit is allowed, but this leads to the extreme dilaton black holes 
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and not to regular particle-like solutions. For such extreme dilaton black holes 
the metric-dilaton relation Na 2 = e 2 ( < t>~ < t>°°)h is restored. At the origin these 
configurations approach the extreme Abelian solutions with the same value of 
charge [138], while showing the usual nodal structure for r > 1. 

The SU(2) EYM(D) black holes admit also axially symmetric generalizations 
[200,202]. By far, these give the only known explicit example of static and 
non-spherically symmetric black holes with non-degenerate horizon. The so- 
lutions can be obtained by extending the spherically symmetric black hole 
configurations to the higher values of the winding number u, similarly to the 
procedure in the regular case. The solutions are asymptotically flat and have 
a regular event horizon, but for v > 1 they are not spherically symmetric and 
the event horizon is not a sphere but a prolate ellipsoid. Specifically, express- 
ing coordinates p, z used in the axial line element (2.21) in terms of r, d as 
z = rcos-$, p = rsin-$, the event horizon is defined by the relation r = r^, 
where (?oo( r fe) — 0. The horizon hence has the S 2 topolgy, but geometrically 
it is not a sphere, since its circumference along the equator turns out to be 
different from that along a meridian. The surfaces of constant energy density 
are either ellipsoidal or torus-like. It is interesting that the energy density is 
not constant at the horizon and depends on ■&. However, the surface gravity 
K sg , where 

K 2 ag = -\g 0o 9 ik d i g Qo d k g Q o, (4.28) 

is constant at the horizon. If a dilaton is included, then the relation (4.26) 
between mass, dilaton charge and the event horizon parameters holds also for 
v > 1. The solutions are characterized by (r/j,z/, n). Here v and n are the 
winding number and the node number of the amplitude w in (2.22), respec- 
tively. The deviation from spherical symmetry increases with growing v, but 
decreases for large n. As n — > oo the solutions tend to the Abelian spherically 
symmetric black hole solutions. In the limit of shrinking horizon, — > 0, 
the solutions reduce to the regular axially symmetric EYM(D) solitons. More 
details and interesting pictures can be found in [202]. 



5 Stability analysis of EYM solutions 

The stability issue for the BK solitons and EYM black holes has been exten- 
sively studied [297,296,343,342,325,323,321,220,70,71,73], both perturbatively 
and at the non-linear level. A novel feature arising in this analysis is the fact 
that, unlike the situation in the (electro)-vacuum case, a non-trivial temporal 
dynamics for EYM fields exists already in the spherically symmetric sector. 
The Birkhoff theorem thus does not apply [19,70], and the physical reason 
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for this can be deduced from the "spin from isospin" phenomenon [178,158]. 
The spin and isospin of the gauge field can combine together giving zero value 
for the sum, in which case the vector YM field effectively behaves as a scalar 
field. As a result, spherical waves exist in the theory. It turns out then that 
all known regular and black hole solutions in the EYM theory with the gauge 
group SU(2) are unstable with respect to small spherically symmetric pertur- 
bations [325]. This is true also for the solutions for higher gauge groups [71,73]. 
The numerical analysis indicates that the non-linear instability growth results 
in the complete dissipation of the initial equilibrium configuration [343]. The 
energy is then partially carried away to infinity by spherical waves and par- 
tially collapses to the center. 

It turns out that the n = 1 BK solution can play the role of the intermediate 
attractor in the gravitational collapse of the YM field [88]. Near the boundary 
in the initial data space between data which form black holes and data which 
do not there is a region characterized by a mass gap. The minimal black hole 
mass for data from this region is equal to the mass of the n = 1 BK solution. 
This is reminiscent of a first order phase transition. There is also another 
region in the parameter space characterized by the echoing phenomenon and 
the scaling behavior for the black hole mass, which can now be arbitrarily 
small [88]. This is explained by the existence of another intermediate attractor 
in the problem. The corresponding solution to the EYM equations is time- 
dependent. It describes the situation when there is a constant ingoing flux 
of the YM radiation which comes from infinity, completely reflects from the 
origin, and then goes back to infinity. The system balances just on the verge of 
collapse but it does not actually collapse due to fine tuning of the parameters. 
The approximate description of this solution was found in [151]. 

In what follows we shall consider only small, spherically symmetric perturba- 
tions for the BK solitons and non-Abelian black holes. The complete set of 
perturbations splits into two parity groups. The even-parity sector is obtained 
via perturbing the non- vanishing gauge and metric amplitudes in Eqs.(2.49) 
and (2.50): 

w — > w + 5w(t, r), m — > m + 5m(t,r), a — > a + 5a(t,r). (5.1) 

The odd-parity group consists of perturbations of those amplitudes in the Wit- 
ten ansatz (2.13) that vanish for static, purely magnetic background solutions. 
These are ao, a r , and Imw: 

a — > 5a (t,r), a r — > 5a r (t,r), Imw — > 5(lmw). (5.2) 

The parity of the modes is determined by Eq.(2.16). Since the background 
solutions are parity-invariant, the two groups of perturbations decouple at the 



56 



linearized level. 



5. 1 Even-parity perturbations 



Substituting (5.1) into the EYM equations (2.44)-(2.48) and linearizing with 
respect to 5w, 5m, and 5a gives the perturbation equations. It turns out that 
the gravitational amplitudes 5m and 5a can be expressed in terms of 5w, which 
leads to a single Schrodinger-like equations for 5w. The procedure is as follows. 
Linearizing Eq.(2.47) one arrives at 5m = 2Nw'5w with N — 1 — 2m/r, whose 
solution with the appropriate boundary conditions is 

5m = 2Nw'5w. (5.3) 

Similarly, linearizing Eq.(2.46) one obtains 

5 (a' /a) = 4w'5w'/r. (5.4) 



Finally, the Yang-Mills equation (2.45) containt after the linearization 5w, 5m 
and 5 (a' /a), where the last two terms can be expressed using (5.3) and (5.4). 
As a result, one arrives at [297] 

(J1 +a ^ N ^_: 2(^1 =^//. (0.5) 



dp 




where 5w(t,r) = exp(iut) i](p), and the functions w, N, and a refer to the 
background solution under consideration. The "tortoise" radial coordinate p 
is defined by 

d P 1 (K R\ 

Tr = ^N> (5 - 6) 



where p G [0, oo) for solitons and p G (— oo, oo) for black holes. It is worth 
noting that Eq.(5.5) can also be obtained via computing the second variation 
of the mass functional M(l,7r/2) in Eq.(3.23) with respect to w [318,316,68]. 
The first variation vanishes on shell, while the second one is 

5 2 M = [ 5w(p) \-^- + U ) 6w(p) alp, (5.7) 



where the potential U coincides with the one in Eq.(5.5). 
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Eq.(5.5) determines a regular Schrodinger eigenvalue problem on a line or 
semi-line. The potential is bounded for black holes, while for solitons it can 
be represented as a sum of a bounded piece and the p-wave centrifugal term 
2/p 2 (see Fig. 11). The numerical analysis reveals the existence of n bound 
state solutions to Eq.(5.5) with u? < for the n-th BK or EYM black hole 
background (with any r^) [297,296]. This implies that the solutions are unsta- 
ble. The non-linear instability growth for the BK solitons has been analyzed 
numerically [343,342]. It has been found that a part of the initial static con- 
figuration collapses forming a small Schwarzschild black hole, while the rest 
radiates away to infinity. It is worth noting that for the EYM black holes with 
n — 1 the negative eigenmode gives rise to the gauge field tensor 5F that is 
unbounded at the horizon [34]. However, it turns out that this mode is nev- 
ertheless physically acceptable, because the divergence can be suppressed by 
making wave packets with real frequency eigenmodes [41,329]. 



5.2 Odd-parity modes 

The odd parity perturbation sector also contains negative modes, whose exis- 
tence can be established [137,47,71,325,73] without even resorting to numeri- 
cal analysis. The number of such modes can also be determined analytically 
[323,321]. For the BK solitons the existence of the odd-parity negative modes 
directly relates to the sphaleron interpretation of the solutions. This is due 
to the fact that the Chern-Simons number of the gauge field changes for pul- 
sations of the type (5.2). The even-parity perturbations, on the other hand, 
preserve the Chern-Simons number. In order to see this it is instructive to 
compute the *FF invariant for the spherically symmetric gauge field (2.13): 



where p\ = Rew and P2 = Imw. This vanishes if ao = a r = Imw = 0. 

The perturbation equations in the odd-parity sector are obtained by linearizing 
the Yang-Mills equations (2. 44), (2. 45) with respect to 5a , Sa r and 5(lmw), 
which gives a system of three linear equations. For time-dependent perturba- 
tions one can impose the temporal gauge condition 5ao = 0. As a result, the 
two independent perturbation equations can be represented in the Schrodinger 
form 



— ti* F /JU/ F'»'=(a 1 (\w 



1) +P2P[ -P1P2)' 




(5.8) 



(5.9) 
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3 



log(p) 



Fig. 11. Potentials U cvcn and ?7 dd in the Schrodinger eigenvalue problems (5.5) and 
(5.14), respectively. 

Here, introducing p = —id /dp and 7* = \ / 2Na/r, the Hamiltonian H and the 
wavefunction \1> are given by 



~ I w 2 ^ 2 ^w' -^wip N 

H = 

K 7*w' p + ip-f*w p 2 + \^l{w 2 - 1) ) 




(5.10) 



with 5a r (t,r) = exp(iut)y2/r 2 N((p) and 5(lmw(t,r)) = exp(iut) £(p) . The 
radial variable p is defined by Eq.(5.6). The third perturbation equation, which 
is the linearized Gauss constraint (Eq.(2.44) with f3 — 0), is a differential 
consequence of (5.9). 

The Hamiltonian in (5.10) is formally self-adjoint and real, which ensures the 
existence of self-adjoint extensions. In order to show that the spectrum has a 
negative part, it suffices to find a function \1> such that 

- 2 W- i 7 ^<0. (5.11) 

Here the scalar product = / ^<&dp. For solitons the function \I/ with 

the required properties is specified by [47] 

C = </7*, Z = w 2 -l. (5.12) 



Note that this can be obtained by differentiating the vacuum-to- vacuum in- 
terpolating family (3.16) with respect to the parameter A at the sphaleron 
position, A = n/2. The direct computation gives 

(mm = -j [w' 2 + 1(^ 2 - 1) 2 ) dp < 0, (5.13) 
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while the norm is finite. An expression similar to (5.12) can be found 

also in the black hole case [325]. As a result, all BK solitons and EYM black 
holes are unstable. 

The above variational argument can be generalized to arbitrary gauge groups. 
The conclusion is: all static, spherically symmetric and purely magnetic regular 
or black hole solutions to the EYM equations for arbitrary gauge groups are 
unstable [71,73]. It is interesting that this conclusion can be reached without 
explicit knowledge of the possible equilibrium solutions. 

The actual number of the odd-parity negative modes for the n-th BK or EYM 
black hole background is n. In order to see this one makes use of the fact 
[323,321] that for to ^ the eigenvalue problem in (5.9), (5.10) can be mapped 
to the equivalent one-channel Schrodinger problem 



+ y ( 3u]2 - x ) + 2 \ 4 = ( 5 - 14 ) 



Here 



V>= (L + ^+tfS^JL , (5.15) 
\dp w I 



and Z is a solution of an auxiliary nonlinear differential equation 

Z'^w'Z 2 -^. (5.16) 



The direct verification shows that for u = Eq.(5.14) has a solution 

tpo = wexp I J w 2 Zdp J . (5-17) 



One can show that Z can be chosen in such a way that the potential in (5.14) 
is everywhere bounded, both for solitons and black holes (see Fig. 11), while ipQ 
is normalizable. Due to the factor w, the number of nodes of ipo is n. The well- 
known theorem of quantum mechanics guarantees then that the eigenvalue 
problem (5.14) has exactly n negative energy eigenstatesP ]. 

In order to actually find the bound states one solves numerically Eqs.(5.9), 
(5.10) or Eq.(5.14) [220,321]. For the n = 1 EYM soliton and black holes 

io p or m ulti-channel Schrodinger problems bound states can be counted by using 
the nodal theorem recently proven in [8]. 
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Table 2 

The eigenvalues of the negative modes for the n = 1 EYM solutions. 



rh 


lo 2 (even-parity) 


(odd-parity) 


(soliton case) 


-0.0524 


-0.0619 


0.1 


-0.0502 


-0.0618 


0.5 


-0.0402 


-0.0600 


1 


-0.0268 


-0.0492 


5 


-0.0021 


-0.0029 


10 


-0.0005 


-0.0007 



the bound state energies for both values of parity are presented in Table 
2. Summarizing, for any n and for all values of the event horizon radius 
the BK solitons and non-Abelian black holes are unstable with respect to 
small spherically symmetric perturbations. The number of instabilities is 2n, 
of which n belong to the odd-parity sector and n are parity-even. 



6 Slowly rotating solutions 

The difference between the non-Abelian EYM theory and its Abelian coun- 
terpart is emphasized further still when one studies stationary generalizations 
of the static BK solitons and EYM black holes [327,69,66]. It turns out that 
rotating EYM black holes acquire an electric charge, and not just a dipole 
correction to the background gauge field which one would normally expect in 
the Abelian theory. As a result, spinning up the neutral static solutions makes 
them charged up. The family of stationary EYM black holes includes in ad- 
dition quite peculiar solutions which are non-static and have non-vanishing 
angular velocity of the horizon, but the total angular momentum measured 
at infinity is zero. This shows that the electrovacuum staticity theorem [78] 
does not generalize in a straightforward way to the EYM theory. The regular 
BK solutions admit charged, stationary generalizations too. This is even more 
unusual, since solitons in other field-theoretical models, such as the t'Hooft- 
Polyakov monopoles, say, cannot (slowly) rotate [66,171]. 

Studying stationary EYM fields presents certain difficulties. Although solu- 
tions for static, axially symmetric deformations of the BK solitons and EYM 
black holes have been obtained numerically [199-202], the problem becomes 
much more involved in the stationary, non-static case [126]. In particular, one 
realizes then that the EYM equations do not imply the Frobenius integrability 
conditions for the Killing vectors [165]. The standard Papapetrou ansatz can 
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therefore be too narrow. This shows that the Abelian circularity theorem also 
does not generalize to EYM systems in a straightforward way. 

Although solutions to the full stationary problem are still lacking, the per- 
turbative analysis based on the assumption of linearization stability has been 
carried out [327,69,66]. The basic idea is as follows. Suppose that there is a 
one-parameter family of stationary black hole or regular solutions of the EYM 
equations, approaching the static solutions for angular momentum J = 0. 
The tangent to this family at J = satisfies the linearized EYM equations. 
Conversely, it is reasonable to expect that for a well-behaved solution of the 
linearized equations around the static configurations there exists an exact one- 
parameter family of stationary solutions. Accordingly, the problem reduces to 
studying the linear rotational excitations for the static BK solitons and EYM 
black holes. 

6.1 Perturbation equations 

Consider small perturbations around the BK solitons or EYM black holes: 

9i*v -»■ 9nu + hfju,, A u -> A v + Vv, (6.1) 

where [g^^AJ) refer to the background configuration, and it is convenient to 
express the gauge field in the string gauge (2.17): 

A = w (T 2 dd — Ti sin dip) + T 3 cos dip. (6.2) 

The perturbation equations are obtained by linearizing the EYM field equa- 
tions (2.4)-(2.6) 

— VffV/i^ — 2 R^ph^ + R°h ai/ + Ryh ali = 4 ST^, 
- D a D a ip v + R^ a - 2[F ua: VI + h aP D a F Pu + F a N a hp v = 0. (6.3) 

Here D v = V v + [A v , ■ ], and the gauge conditions 

V„K = K = D„r = (6.4) 

are imposed. 5T^ U in Eqs.(6.3) is obtained by varying the energy-momentum 
tensor 

= X - tr (F^ a F vP g a " - \g, v F aP F pfJ g a Pg^ (6.5) 
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with respect to the metric and the gauge field. 



6.1.1 General mode decomposition 

In order to identify the most general rotational degrees of freedom, one de- 
termines those amplitudes in the partial wave decomposition which can give 
a non-vanishing contribution to the ADM flux integral for the total angular 
momentum 

s 2 



In order to carry out the partial wave decomposition for (h^,tp u ), it is con- 
venient to introduce the complex 1-form basis 9 a : 

e° = aVNdt, e 1 = -^=, e 2 = -^(dd-i sine dp), e 3 = (e 2 y, (6.7) 

y N v 2 



the non- vanishing components of the tetrad metric r] a ^ = (8 a , 0") being 77°° = 
— r/ 11 = —i] 23 — 1. In addition, one introduces the new Lie-algebra basis: 
Li = T x + iT 2 , L 2 = T x — iT 2 , L 3 = T 3 . The perturbations then expand as 

V = ^H a/3 , ^ = L a *£0» (6.8) 



and the complete separation of variables in the perturbation equations (6.3) 
is achieved by making the following ansatz: 

H a/3 = exp(iwi) H af3 (r) s Y jm (#, <p), V a a = exp(iwi) $£(r) 3 Y jm (&, <p), (6.9) 



where s Yj m { r d ) p) are the spin-weighted spherical harmonics [148]. Here the 
quantum numbers j and m are the same for all amplitudes. The spin weight 
s can be different for different harmonics and is given by s = s a p for a a /3, 
where 

s 00 — s 01 — s ll — s 23 — 0, S02 — s 12 — — s 03 — ~~ s 13 = 1? s 22 = — s 33 = 2, 

and by s = with 

s = S l = S 3 = 0' s = S l = S 2 = ~ s = ~ s l = ~ s 3 = 1) S 2 = ~ s 3 = ^ 

for the YM perturbation amplitudes VP® . As a result, for given u and j > in 
(6.9), Eqs.(6.3) reduce to a system of radial equations for H a p(r) and 
Due to the spherical symmetry of the background fields, the quantum number 
m does not enter the radial equations. 



63 



6.1.2 Rotational modes 



Inserting (6.9) into (6.6) the integral is computed by integrating over a two- 
sphere at finite radius r and then taking the limit r — > oo. The angular 
dependence of the integrand implies that the integral vanishes for any r unless 
j — 1, and that only the h 0tp perturbation component can give a non- vanishing 
contribution. Now, the transformation behavior of the angular momentum 
under space and time reflections (P,T) implies that only those perturbation 
amplitudes are relevant which are even under P and odd under T. For j = 1 
these appear only in ho v and in two isotopic component of ipo. They decouple 
from the remaining modes because the background solutions are P and T 
symmetric. This leads finally to the following most general ansatz (up to global 
coordinate rotations) for the stationary rotational modes: 

h = 2S{r) sin 2 6 dtd V , $ = ( T l ^- sin Q + T 3 ^ cos 6 j dt . (6.10) 



The gauge conditions (6.4) for this ansatz are fulfilled identically and the 
perturbation equations (6.3) reduce to the coupled system for the radial am- 
plitudes S, x, and m 



+ 4iVrV (j^V 4 ^ 2 ^ X) (w X -v) = 0, 
-r 2 Na (^j + (l + w 2 - 2w' 2 N) X 

-2w V -rN (w'S)' + (^2Nw' 3 + S = 0, 

-j^Na^pj + 2 (l + w 2 - w' 2 N) v - Aw X + 2(1 ~ — S = 0. (6.11) 

Here w, m, and a refer to the background solutions, N = 1 — 2m/ r. Given a 
solution of these equations, the ADM angular momentum is 

J l = ^}^(^)'. (6.12) 



r— ►oc \6r 



Note that, with a suitable reparameterization of the variables S, x-> an( i V-> 
Eqs.(6.11) can be represented in a Schrodinger form with a manifestly self- 
adjoint Hamiltonian [66,69]. 



64 



6.2 Solutions 



For the Schwarzschild background, N — 1 — 2M/r, a = w = 1, Eqs.(6.11) 
admit the solution x — V — 0, S = —UMjr. This is recognized as the linear 
rotational excitation of the Schwarzschild metric. One can see that this mode 
is bounded everywhere outside the horizon. However, since it does not vanish 
at the horizon, it is not normalizable. Accordingly, when looking for solutions 
to Eqs.(6.11) for the non-Abelian backgrounds, one considers all modes that 
are bounded and regular, but not necessarily normalizable. The strategy then 
is to find all well-behaved local solutions in the vicinity of the origin or horizon 
and at infinity, and then match them in the intermediate region. 

In the far field region the most general solution to (6.11) that gives rise to 
regular perturbations in (6.10) reads 



x = Vr + Q + 0(r~ 1 ), rj = -Vr - Q + O^ 1 ). (6.13) 

Here a, M are the parameters entering the asymptotic expansion (3.6) for the 
background w, m. This solution contains four integration constants: 3, V, Q, 
while the fourth one, C4, is contained in the higher order terms in (6.13). 

In the near zone, in the soliton case, the most general local solution to (6.11) 
that is regular at the origin is 

S = Cl r 2 + 0(r 4 ), x = c 2 r + 0(r 3 ), rj = c 2 r + 0(r 3 ). (6.14) 



This contains three independent parameters: Ci, c 2 , while the third one, C3, 
resides in the higher order terms in (6.14). Since the coefficients in (6.11) are 
continuous and regular for < r < 00, the local solutions (6.13) and (6.14) in 
the vicinity of r = and r = 00 admit extensions to the semi-open intervals 
[0,oo) and (0,oo], respectively. The total solution space is six dimensional, 
while the local solutions specify four and three dimensional subspaces, whose 
intersection is (at least) one dimensional. In other words, the matching condi- 
tions for S, x, an d rj give six linear algebraic equations for the seven coefficients 
J, Q, V, c\ . . . C4, such that the matching is always possible. For example, for 
the n = 1 BK solution the numerical matching gives 

Q/J = 0.960, Vj 3= -0.594, (6.15) 



and similarly for c\,...Ci. All BK soliton solutions therefore admit a one 
parameter family of stationary excitations parameterized by 3. 
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For black holes the space of solutions that are regular in the vicinity of the 
horizon is four dimensional 



S = c r h + ci x + 0(x 2 ), 

X = -c w(r h ) + c 2 x + 0(x 2 ), 7] = -c + c 3 x + 0(x 2 ), (6.16) 

where x — r — Th and Co, ci, C2, and C3 are four independent integration 
constants. As a result, all EYM black holes admit a two parameter family of 
stationary excitation. 

In order to clarify the meaning of the parameters in the solutions, one uses 
Eq.(6.12), which shows that the constant J is the ADM angular momentum. 
Next, after the gauge transformation with U= exp(i(7r — i?)T 2 ) the gauge field 
perturbation in the asymptotic region becomes 

i> = T 3 (v + Q) dt + 0(r' 2 ). (6.17) 



This suggests that Q should be identified with the electric charge. The con- 
stant V determines the asymptotic value of the temporal component of the 
gauge field, Aq(oo). In the Abelian theory it would be possible to gauge V 
away. In the non- Abelian theory, however, such a gauge transformation would 
render the whole configuration time- dependent. The physical significance of 
V became clear already in the course of the matching procedure. In addition, 
due to the coupling to the background gauge field, it enters the asymptotic 
expansion of the field strength, SF, whose non-vanishing components in the 
same gauge as the one used in (6.17) read 

5F 0r = ^T 3 + O(r- 3 ), 5F # = ^T 1 + O(r- 2 ), 

§F 0v = -— T 2 sintfcostf + 0(r- 2 ). (6.18) 
r 

It is obvious that Q can be expressed in terms of the flux of the electric field 
over a two-sphere at spatial infinity. The constant V, on the other hand, gives 
no contribution to the flux, since the corresponding piece of electric field in 
(6.18) is tangent to the sphere. 

Summarizing, all BK soliton solutions admit slowly rotating excitations with 
continuous angular momentum J and electric charge Q proportional to J. 
Note that solitons in other field theoretical models, such as the t'Hooft- 
Polyakov monopoles or boson stars, generically do not admit slow-rotating 
states [327,69,66,171]. This might be due to the fact that a soliton, being a 
solution to a non-linear boundary value problem, can exist only for discrete 
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values of the parameters. The angular momentum should therefore be quan- 
tized. In this sense the situation in the EYM theory is exceptional. This is 
presumably because both gravitational and gauge fields are massless, which 
manifests in the slow (polynomial) decay of the fields at infinity. 

The slow rotational excitations of the EYM black holes can be parameter- 
ized by any two of the three parameters J, Q, or V. Accordingly, there are 
three distinguished branches of solutions. First, for solutions with V — the 
charge is proportional to the angular momentum [327]. Secondly, there is an 
uncharged branch with Q = 0. Finally, there are solutions with J = 0, which 
are non-static, as can be seen already from (6.13). In view of this the Abelian 
staticity theorem asserting that stationary black holes with J = must be 
static does not apply. The non-Abelian version of this theorem [300] states 
that a stationary EYM black hole is static and the electric field vanishes, 
provided that the following condition holds 

Q H J-VQ = 0, (6.19) 

where fin is the angular velocity of the horizon. Notice that this agrees with 
the results described above. Indeed, (6.19) implies that stationary black holes 
with J = must be static if only VQ = 0, which is not the case. On the 
other hand, for non-static solutions the condition (6.19) should be violated, 
which, in particular, shows that rotating solitons must be charged up, since 
with n H = one has to have VQ ^ 0. 



7 Self-gravitating lumps 

In this chapter we shall discuss gravitating solitons and black holes in the theo- 
ries admitting particle-like solutions in the flat spacetime limit. Historically, 
the first investigation of such systems was motivated by a wish to understand 
the structure of very heavy magnetic monopoles [312]. It was found that the 
gravitating generalizations for flat space monopoles exist at least for small 
values of Newton's constant. More systematic investigations of the problem 
were undertaken after the discovery of the BK solutions. These have revealed 
a number of the typical features. First, it has been found that, for a large 
class of non-linear matter models, the flat space solitons can be generalized 
to curved spacetime, provided that the dimensionless gravitational coupling 
constant k is small. This is not very surprising, since it is intuitively clear 
that for small values of k perturbation theory applies. However, and this is far 
less obvious, it turns out that under certain conditions gravity can be treated 
perturbatively even for black holes, provided that the event horizon radius 
r h is small. As a result, the gravitating lumps can be further generalized by 
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replacing the regular center by a black hole with a small radius tv Below 
we shall review the corresponding argument based on the implicit function 
theorem [188] . This is so simple that the existence of hairy black holes might 
actually have been foreseen many years ago. 

Secondly, in the non-linear regime, when gravity is not weak, the fundamental 
lumps and black holes admit a discrete spectrum of gravitational excitations. 
The excited solutions become infinitely heavy for k — > and, remarkably, re- 
duce then to the rescaled BK solutions. The excitations thus can be thought 
of as gravitating solitons with small BK particles or EYM black holes in the 
center. Finally, the fundamental and excited lumps and black holes cannot 
exist beyond certain maximal values of k and r^. The existence of a bound for 
k is easy to understand. As k is proportional to the ratio of the gravitational 
radius of the object to its typical size, it cannot be too large, since otherwise 
the system becomes unstable with respect to the gravitational collapse. The 
solutions, however, do not collapse as k approaches the critical value. Instead 
they either become gravitationally closed or coalesce with the excited solu- 
tions. The existence of the bound for r\ is quite interesting. It seems as if a 
small black hole could not swallow up a soliton which is larger than the black 
hole itself. As a result, a hairy black hole appears. However, a big black hole 
can completely absorb all non-trivial hair. 



7.1 Event horizons inside classical lumps 



One can use the implicit function theorem to argue that, within a large class of 
non-linear models, and for small values of the gravitational coupling constant k 
and the event horizon radius r^, the flat space solitons admit both regular and 
black hole gravitating generalizations. This comes about as follows [312,188]. 
Consider a static, spherically symmetric gravitating system with the metric 
(2.50) and the stress tensor Tff = diag (p, — p r , —p e , —pe). With the notation 
of Sec. 2, the non-trivial Einstein equations G^ u = 2kT^ u are 

m' = r 2 p, a' = k — (p + p r ) a, (7.1) 



with N = 1 — 2/tm/r. This has to be supplemented with the condition V^T^ 1 = 
0, which can be rewritten with the use of (7.1) as the Oppenheimer-Volkoff 
equation 

Pr = - K r 2 N " (P + Pr) + -(P0-Pr)- (7.2) 



Suppose that there is an equilibrium flat space configuration, such that equa- 
tions (7.1) and (7.2) are fulfilled for k — 0. In most cases one can expect 
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that this configuration will also survive if k is non-zero and small. Indeed, for 
k < 1 the additional terms in the equations are small and gravity can be 
treated perturbatively. One will have then N m a ~ 1. At the same time, it is 
less clear whether the corresponding black hole generalization will also exist. 
For black holes N vanishes at r h implying that the deviation from the flat 
space value N = 1 is non-small. This also implies that one cannot use per- 
turbation theory, since the coefficient 1/N in (7.1), (7.2) diverges at r = tv 
However, if the matter model is such that (p + p r ) ~ iV at the horizon then 
the blowing up of 1/N will be canceled, in which case gravity can be treated 
perturbatively even for black holes [188]. 

It has already been mentioned that for all matter models considered in Sec. 2. 3 
the reduced 2D matter Lagrangian in the static, purely magnetic case can be 
represented in the form 

L m = -(NK + U), (7.3) 



where K and U depend on r and on the matter field variables, but not on the 
gravitational variables iV and a. For this Lagrangian one has 

p = (NK + U)/r\ p T = (NK - U)/r 2 , (7.4) 



such that (p+Pr) = 2NK vanishes at the horizon. As a result, given a flat space 
soliton solution, one can expect that it will have gravitating generalizations, 
both in the regular and black hole cases. The argument is as follows. Inserting 
(7.4) into (7.1) the Einstein equations for m and a become linear, and hence 
can be integrated in quadratures for given K and U . For example, one will 
have 

(j = exp {-2k J K ~^j ( 7 - 5 ) 



and similarly for m. Therefore, in order to solve the gravitating problem, it 
remains to determine the matter variables entering K and U. Let us denote 
these variables collectively by <f>. The equations of motion for can be obtained 
by varying the ADM mass functional [312]. The ADM mass is M = ra(oo) = 
m(oo)er(oo). Using (ma)' — (K + U)a one obtains 

oo / oo \ 

M[0, k, r h \ = T ± a(r h ) + J dr(K + U) exp I -2k J K . (7.6) 

rh \ r / 



This expression takes gravity into account but depends only on the matter 
variables; notice that for r h = it was used in Eq.(3.23). Varying M with 
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respect to gives the matter equation of motion with all gravitational degrees 
of freedom expressed using formulas like (7.5). The problem therefore reduces 
to studying solutions of 

F[<j>,K,r h ] = ^ = 0. (7.7) 



By assumption, there exists a flat space solution <po, such that F[cf>o, 0, 0] = 0. 
Suppose that the operator SF/S<p is invertible for this solution. This will be the 
case, for example, if the solution is a local minimum of energy, 5F[cf) , 0, 0] > 0. 
Then the implicit function theorem ensures that for k and sufficiently close 
to zero, there exists 0(/c, r h ) satisfying (7.7) such that 0(0,0) = O [188]. 

Summarizing, if a flat space soliton solution fulfills the conditions specified 
above then it will survive also in the weak gravity case, and in addition 
will be able to contain a small black hole inside. This happens, for exam- 
ple, for monopoles, sphalerons, and Skyrmions considered below. However, 
not all matter models are of the type (7.3). For example, systems with electric 
fields are not of this type and hence are not covered by the above arguments. 
Another example is provided by boson stars. For the BK particles there are no 
flat space counterparts, and the above arguments do not directly apply. How- 
ever, it follows that for k ^ the existence of regular solutions to F[<f>, k,0] = 
implies the existence of corresponding black hole solutions to F[<f>, K,r h ] = 0. 
As a result, the existence of the BK solitons implies the existence of the EYM 
black holes, at least for small r h . 



7.2 Gravitating monopoles 



The regular monopole solutions were discovered by t'Hooft and Polyakov in 
the YMH model (2.32) with the triplet Higgs field [305,263]. For monopoles 
in flat space there is the absolute lower bound for energy determined by the 
topology of the Higgs field [44]: 

Attv 

E>—\v\. (7.8) 
g 

Here v is the winding number of the Higgs field, which coincides with the 
magnetic charge. Solutions saturating this bound are topologically stable. For 
the hedgehog ansatz in (2.33) one has v — 1, and the corresponding EYMH 
equations following from (2.34) in the static, purely magnetic case read 



m' = Nw' 2 + ^- + ^N(j)' 2 + w 2 (j) 2 + r 2 V(0), (7.9) 
2r 2 2 
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(AW)' 




(7.10) 



(7.11) 



(r 2 Na(f)'y 



a (2u> 2 + rV(0)) 



(7.12) 



with 



N =1- 



2 nrn 



^(0 2 -i) 2 - 



(7.13) 



r 



Here k = 4nGv 2 and the length scale is L = 1/gv = 1/M W . Note that ex can be 
eliminated from the equations in the usual way. For a given solution to these 
equations the total energy is the ADM mass M = ra(oo). The dimensionful 
energy is (4irv / g) m(oo). The boundary conditions for finite energy solutions 
are 



w = l-br 2 + 0(r A ), <P = cr + 0(r 3 ), m = 0(r 3 ), (7.14) 
iu = S/iexp(-r), 0= l + C/ 2 exp(-er), m = M + 0(r-- 1 ) (7.15) 



at the origin and at infinity, respectively. Here b,c,B, C, and M are integration 
constants, and the leading behaviour of the functions fi(r) and /2(f) for r — > 
00 is power law [89,192]. 

7.2.1 Flat space solutions 

The flat space limit corresponds to k — * 0, in which case the YMH equations 
(7.11) and (7.12) decouple, while the Einstein equations (7.9) and (7.10) give 
N = a = 1. The mass of the solutions, m(oo), is obtained from (7.9), which 
gives after simple rearrangements 



Using (7.14) and (7.15) the integral of the last term in the integrand is equal 
to one, while the first three terms are positive definite. This shows that 



M 




(7.16) 



M > 1, 



(7.17) 
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in agreement with the bound in (7.8) for v = 1. In fact, (7.8) is obtained in 
a similar way by picking a total derivative in the 3D mass functional. Now, if 
e = and 

w' + w<p = 0, r 2 0' + w 2 -l = (7.18) 



then the first three terms in the integrand in (7.16) vanish and the bound 
M = 1 is saturated. The solution of the Bogomol'nyi equations (7.18) describes 
the BPS monopole [265]: 

r 1 
w = — — -, = coth(r) - -. (7.19) 
sinh(r) r 



By construction, this solution is stable. It has only one non-vanishing compo- 
nent of the stress tensor, T °, and hence the condition in (3.1) is fulfilled. 

For e > the bound in (7.17) is not attained and the first order equations 
cannot be used. One has to solve then the second order YMH equations (7.11) 
and (7.12). Solutions with the boundary conditions (7.14), (7.15) can be ob- 
tained numerically [45,193]. They describe the t'Hooft-Polyakov monopoles. 
For any given value of e there is one solution, the existence of which was 
shown in [311]. The uniqueness of the e = solution was established in [233]. 
The mass M increases with e such that 1 < M(e) < M(oo) = 1.96. In the 
limit e — > oo the Higgs field becomes "frozen" and there remains only the YM 
equation (7.11), where one sets 0=1 [193]. The t'Hooft-Polyakov monopoles 
are linearly stable [12]. For any e the monopole solutions can be generalized 
to include an electric charge [183]. 



7.2.2 Regular gravitating monopoles 

For k^O the energy condition (7.17) no longer holds due to the gravitational 
binding. In fact, for a given k, taking the boundary conditions into account, 
the mass functional M[0, k, 0] in (7.6) with 

K = w" + ^<P'\ U=^^+w^ + ^r 2 (0 2 -!) 2 , (7.20) 



is still bounded from below by a non-zero value [312]. However, it is unclear 
whether the M min (/t) = inf M[0, k, 0] is reached in the set of solutions of the 
EYMH equations. For example, it turns out that if k is large, the fields that 
minimize M[0, k, 0] do not even belong to the space of differentiable functions 
[222]. At the same time, the argument of the preceding section suggests that 
the gravitating monopoles exist at least for small values of k. This shows that 
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Fig. 12. On the left: Amplitudes w n , <f> n , and N n for the fundamental, n = 0, 
monopole solution (solid lines) and for its first, n = 1, excitation (dashed lines) for 
re = 0.6 and e = 0. On the right: the ADM mass versus re for the n = and n > 1 
regular monopole solutions with e = 0; and the mass Mrn = Vv^ °f the extreme 
RN solution. All curves for n > 1 terminate at the point B with re = 3/4, while 
staying very close to each other and to the RN curve for re < 3/4. For small re the 
mass of the n-th solution is M n / y/k, where M n is the mass of the n-th BK soliton. 
The n = curve terminates at re = 1.97. At the point B the extreme RN solution 
bifurcates with infinitely many non-Abelian EYMH black hole solutions. 



one can study these solutions numerically, however, the issue of their stability 
remains open. 

The numerical integration of Eqs.(7.9)-(7.12) with the boundary conditions 
(7.14)-(7.15) reveals the following picture [258,222,54,5,56,232]. First, for re 
being small enough, there are the self-gravitating monopole solutions. These 
are qualitatively similar to the flat space t'Hooft-Polyakov monopoles and 
reduce to them for re — > 0. The amplitudes w, <p, m, and a are monotone 
functions, while N develops a minimum at some r m ~ 1, N(r m ) > 0. One can 
call r m the monopole radius. The solutions are characterized by re and e. The 
mass, M(re), decreases with growing re due to the gravitational binding (see 
Tab. 3 and Fig. 12). At the same time, the monopole radius does not change 
considerably. As a result, the ratio of the gravitational radius of the monopole 
to its radius r m is proportional to re. For regular solutions re must not therefore 
exceed some critical value, otherwise the system becomes unstable with respect 
to the gravitational collapse. 

In agreement with the last remark, the self-gravitating monopoles exist only 
for a finite range < re < re max (e). One has re max (oo) = 1/2 < re max (e) < 
K max(0) = 1.97 [56]. As re approaches the critical value, N(r m ) tends to zero 
(see Fig. 12), the functions w and reach their asymptotic values already at 
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r = r m , and the proper distance 




(7.21) 



diverges for r — > r m [^]. As a result, the spatial geometry on the hypersurface 
t =const. develops an infinite throat separating the interior region with a 
smooth origin and non-trivial YMH field from the exterior region where w = 

0, = 1, and the metric is extreme RN. The throat is characterized by a 
constant radius, r Q = y/n max . The metric function cr(r), normalized by a(0) = 

1, diverges for r — > r m . The limiting solution thus splits into two independent 
solutions, which is similar to the behaviour of the BK solutions for n — ► oo. 
The interior solution is geodesically complete. 

7.2.3 Gravitational excitations 

For the gravitating monopoles there is an infinite sequence of radial excita- 
tions. These exist for n being small enough, but do not admit the flat space 
limit, since their masses diverge for k — > 0. The boundary conditions (7.14)- 
(7.15) are still fulfilled, but the amplitude w now oscillates (see Fig. 12) and 
its node number (n = 1, 2 . . . ) characterizes the solutions. This similarity with 
the BK case is not accidental. In fact, for k —>■ 0, the excitations correspond 
to the rescaled BK solutions. This comes about as follows. When k = AirGv 2 
tends to zero, this can be understood either as the weak gravity limit, G — > 
with fixed v, or as the limit v — > with fixed G. Notice that in the latter case 
the length scale L = 1/gv diverges. Consider, however, the rescaling 

r — > \fnr, m^m/y/n, <f> — > (j)/ \/~K . (7.22) 

The whole effect of this on the EYMH equations is that N and V{4>) in (7.13) 
are replaced by 

N=l-— , V{<P) = Utf-Kf. (7.23) 
r 4 

As a result, the limit k — > corresponds now to V(<f>) vanishing (since 2 < k), 
and the EYMH system (7.9)-(7.12) reduces to the EYM equations, which 
leads to the appearance of the BK solutions. Restoring the original length 
scale, the excited solutions for small k consist of the very small (Planck size) 
BK particles sitting inside the large (1/M W ) size monopole. The amplitude w 
oscillates n times in the short interval near the origin and then tends to zero 

11 For small e the picture is slightly more complicated; details can be found in [54]. 
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Table 3 

Masses of the fundamental monopole and its first excitation (e = 0). 



n\n 





0.02 


0.1 


0.5 


0.75 


1 


1.5 


1.96 





1 


0.996 


0.985 


0.926 


0.887 


0.854 


0.783 


0.714 


1 


oo 


6.629 


3.108 


1.413 


1.154 









for r — > oo. For k — ■> the mass tends to M n /^/~K, where M n is the mass of 
the n-th BK solution. All excitations have the same value of /t max = ^max(e)- 
One has n max (0) = 3/4 > « max (e) > /t m ax(oo) = 1/2 [56]. For 

^ * K max I C J 

the excited solutions show the same limiting behaviour with a throat as the 
fundamental solution. 

The picture described above remains qualitatively the same for any < e < oo. 
A novel feature arises for e — > oo [54,5,56]. The Higgs field is frozen in this 
case, and the EYMH field equations reduce to those of the gravitating gauged 
non-Abelian 0(3) sigma-model. Eq.(7.12) then should be dropped, and there 
remain three equations (7.9)-(7.11), where one sets = 1. Note that these 
admit the first integral (7.29). The boundary conditions at the origin for w 
and N are given by 

w = l + 0(r 2 ), N = 1 — 2k 2 + 0(r 2 ), (7.24) 



where k G [0, 1/2]. Since N(0) ^ 1, there is the solid angle deficit and a conical 
singularity at the origin. This is due to the fact that the energy density Tq 
contains the term w 2 (f> 2 /r 2 , which blows up, because does not vanish at the 
origin. For r > the solutions for w and N are qualitatively similar to those 
for e < oo. In particular, they are asymptotically flat (N(oo) = 1) and exist 
for arbitrary node number n. The existence of these solutions was rigorously 
established in [56]. 

In this connection it is instructive to mention also the gravitating global 
monopoles [15]. These arise in the coupled Einstein-Higgs model, which corre- 
sponds to the EYMH theory with dynamics of the YM field being suppressed. 
The field equations in the spherically symmetric case can be obtained from 
(7.9)-(7.12) by ignoring the YM equation (7.11) and setting w — 1 in the 
remaining equations. The boundary conditions for and m at r = are still 
given by (7.14) (if only e < oo), and the origin is regular: N(0) = 1. However, 
since — > 1 for r — > oo while w — 1, the term w 2 <p 2 /r 2 in Tq leads to the diver- 
gence of the energy in such a way that iV(oo)=const.< 1, which corresponds 
to the solid angle deficit at infinity. Note that the absence of finite-energy 
solutions in the model agrees with the general arguments in Sec. 3. 
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For e — > oo the Einstein-Higgs theory reduces to the global gravitating 0(3) 
sigma-model. In this case there remain only two non-trivial field equations, 
(7.9)-(7.10), where one sets w = (p = 1. The solution is a = 1, N = 1 — 2k, 
such that there is a constant solid angle deficit. It is worth noting that this 
model admits the Bogomol'nyi bound and non-trivial exact solutions in the 
case of cylindrical symmetry [98]. 



7.2.4 EYMH Black holes 

For black holes, instead of (7.14), one has the boundary conditions at the 
regular horizon, which are given by the expression similar to that in (4.25) with 
w' h , <p' h , and N' h being functions of Wh and 4>h- The possibility of a degenerate 
horizon was discussed in [152,153,104,32,56]. According to the argument of 
the preceding section, one can expect the black hole generalizations of the 
regular monopoles to exist at least for small values of r^. This is confirmed 
by the numerical analysis for the fundamental (n = 0) monopole and for its 
excitations (n > 1). For <C 1 the solutions resemble small black holes 
sitting in the center of the regular lumps, the latter being almost unaffected 
for r > by the presence of the black hole. The solutions exist only for a 
limited region of the (k, r^) space. The domain of k can be subdivided into 
two parts: < k < 3/4 in which case solutions for any n are possible, and 
3/4 < k < K max = 1.96, which leads only to the fundamental solutions. In 
both cases the solutions exist if only Vh is less than a maximal value, fh{K), 
where /^(re) < 7^(3/4) = Vo/2. For r h — > the black hole solutions tend to 
the corresponding regular ones. In the opposite limit, — > f/j(/«), black holes 
with k > 3/4 develop the infinite throat and become gravitationally closed. 
The metric amplitude N in this case, apart from a simple zero at r = r^, 
develops also a double zero at r = y/n > r^. For k < 3/4 the picture is more 
complicated. In this case for some values of k and one finds more than one 
non-Abelian black hole solution. In addition, for r h > yfn and M > l/y/R 
there exist Abelian RN solutions: 

2kM k 

w = 0, 0=1, a = l, N=l + — . (7.25) 



As a result, the no-hair conjecture is violated. When changes, different non- 
Abelian solutions bifurcate with each other or with the RN solution [54] . This 
leads to the existence of the upper bound for r^. The value f^(3/4) = v^3/2 is 
the least upper bound for for all solutions (for any n and e) with k — 3/4. 
All of them tend for — > Vo/2 to the same limiting configuration, which is 
the extreme RN solution. A more detailed discussion can be found in [54,56]. 



The numbers are given for the case where e = 0. 
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Note that the existence of an upper bound for the black hole radius is quite 
typical for gravitating lumps [256]. One can think that a small black hole is 
unable to swallow up a monopole whose radius is much larger than r^, however 
a bigger black hole can do this, such that all non-Abelian structures disappear 
inside the horizon. This phenomenon does not exist, however, in the pure 
EYM theory, where there is no energy scale other than Planck's mass. The 
non-Abelian EYM black holes can be arbitrarily large. It is worth noting also 
that the fundamental and excited gravitating monopoles and black holes can 
be generalized to include an electric charge [63]. 

7.2.5 Deformed EYMH black holes 

The fundamental gravitating monopoles and black holes are stable, while the 
excited solutions are unstable with respect to small spherically symmetric per- 
turbations. This was shown in [173] with the use of the Jackoby criterion. The 
RN solution (7.25), which is stable in the Abelian theory, becomes unstable 
in the non-Abelian case for k < 3/4 [226,41,221,54,5,56]. This comes about as 
follows. Perturbing the amplitudes w, m, and a of the RN solution (7.25) as 
in Eq.(5.1), while <ft —> <f) + 5<f)(t,r), and linearizing the full system of EYMH 
equations, the linearized YM equation decouples from the rest of the system. 
It is convenient to rescale r and M in (7.25) according to (7.22), which leads 
to the following equation for 5w(t,r) = exp(iut) r](p): 



where N = 1 — 2M/r + 1/r 2 and dp = dr/N. Notice that for k — > 0, which 
now corresponds to the limit where the Higgs field decouples, this equation 
reduces to that in (5.5). As was observed in [56], in the extreme RN case, 
where M = 1, this equation admits infinitely many bound states for k < 3/4, 
because the corresponding potential for p — > — oo behaves like (k — l)/p 2 . For 
k > 3/4 there are no bound states [56]. This corresponds to the fact that for 
k, = 3/4 the extreme RN solution bifurcates with infinitely many non-Abelian 
black hole solutions [56]^. For non-extreme RN solutions the situation is 
more complicated, but the result is the same: the instability arises only for 
small enough values of M and k. 

The instability of the RN solution within the context of the non-Abelian 
EYMH theory has important implications, since it can be viewed as indi- 

13 The exterior region of the extreme RN solution can also coalesce with the corre- 
sponding part of the limiting solutions with the throat. This happens for n > 3/4 
(see Fig. 12). However, since the interior parts are different, this does not lead to the 
change of stability, and the RN solution for k > 3/4 is stable. 




(7.26) 
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cation of the existence of new non-Abelian solutions [223,273,271,272,332]. 
This basic idea is as follows. First of all, notice that the RN solution in (7.25) 
can be generalized to arbitrary integer values of the magnetic charge. For this 
one uses the ansatz (2.24) for higher winding numbers is, while the Higgs field 
is $ = 0T3. The solution is given by 

2kM kis 2 

a = w = 0, 0=1, (7=1, N=l — + — , (7.27) 



and the gauge field strength is F = T3 v sin 1? di? A dip. For any given is this RN 
solution can be either stable or unstable - depending on values of k and M 
[223,272]. Suppose that k and M are chosen such that the solution is close to 
the border of instability and is only "barely" unstable. One can expect then 
that there is a nearby stable non-Abelian solution with lower mass that bifur- 
cates with the RN solution when k and M approach the border. It is plausible 
that this new solution differs only slightly from the RN one, and can be ap- 
proximated with the use of linear perturbation theory [273,271,332]^. Now, 
for v = 1 the static perturbations around the RN solution are spherically sym- 
metric and approximate the non-Abelain EYMH black holes described above. 
However, for is > 1 the solutions of the perturbation equations are no longer 
spherically symmetric^] and exhibit a complicated dependence on spherical 
angles. As a result, assuming that the linear modes indeed approximate some 
new solutions, these new black holes turn out to be not even axially symmet- 
ric. Their shape can be quite complicated, exhibiting discrete, crystal-type 
symmetries (see [273] for interesting pictures). 

It is worth noting that the deformed black holes have actually been described 
within a more general theory, which includes the EYMH model as a special 
case [223,273,271,272,332]. In this theory the third isotopic components of the 
gauge field is regarded as an Abelian vector fields, A^, while the first two 
component are viewed as a complex vector field, W^. The gravitational part 
of the action is standard, while the matter Lagrangian is 



C = - - T^T*" - - H; V H^ + a 

-bd^ + m 2 (0) W*W + l - d^dy - V{4>). (7.28) 

Here = <9 M A - d v A^ and = D^W U - D V W^, while d M „ = i{W*W u - 
W*W,j). The covariant derivative is D^Wy = (<9 M — iA^)W u , and a, b are 

14 The corresponding perturbation equations are not quite the same as in the stan- 
dard perturbation analysis [273]. 

15 It has already been mentioned that all spherically symmetric YM fields with is > 1 
are necessarily embedded Abelian. It follows then that non-Abelian fields for is > 1 
cannot be spherically symmetric. 
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parameters. For a = 1/2, 6 = 4, and m(0) = <fi the theory reduces to the 
triplet EYMH model in the string gauge (2.17) with $ = 0T 3 . For a range 
of the parameters the model admits soliton and black hole solutions, and, in 
particular, the deformed black holes described above. Together with the results 
obtained in the pure EYM theory [200,202], this suggests that the existence of 
deformed static black holes is typical for the gravitating non-Abelian models. 

7.2.6 The Bogomol'nyi bound 

Unfortunately, the flat space Bogomol'nyi equations (7.18) do not generalize 
to curved spacetime. The existence of the Bogomol'nyi bound in flat space has 
a very deep physical reason: the YMH theory for e = can be obtained via 
truncation of a supersymmetric model. On the other hand, the EYMH theory 
is not a truncation of a supergravity model, and no Bogomol'nyi bound exists 
in this case. As a result, a search for Bogomol'nyi equations for k^O leads to 
an almost trivial result: one can find first order equations, but they admit only 
the extreme RN solution [97,14]. For e — > oo, when the Higgs field is frozen, 
the Bogomol'nyi equations are less trivial [124, 53][" re "|: 

k = 1, Nw' = w, r 2 (o-N)' = a{w 2 - 1), rN + w 2 - 1 = ±r. (7.29) 

However, all their solutions, apart from the extreme RN one, have naked 
singularities. In order to obtain a physically meaningful Bogomol'nyi bound 
in curved space [141], one should extend the EYMH theory to make it a part 
of a supergravity model. Such an extension was considered in [142,146]. In 
addition to the EYMH fields it contains a dilaton and an abelian vector field, 
which is typical for toroidally compactified string theory models. The BPS 
solutions in this case can be obtained analytically. 

Sometimes gravitating monopoles are associated with inflation [225,315]. This 
is due to the fact that, since the Higgs field vanishes for r — > 0, there is a false 
vacuum in the monopole core. The results of the numerical analysis of the 
time-dependent problem suggest that monopoles indeed inflate if k exceeds 
the critical value, k > K max [279,278]. 

Other related topics we would like to mention include solitons for the EYMH 
system with the ghost Higgs field [93] , monopole solutions on the cosmological 
background [244], and the Rubakov-Callan effect for the colored black holes 
[132]. 



6 The existence of these equations has no explanation at present. Note that in fiat 
space there are no Bogomol'nyi equations for e — > oo. 
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7. 3 Gravitating YMH sphalerons 

Consider the EYMH model (2.35) with the doublet Higgs field. The spherically 
symmetric ansatz (2.36) contains two independent amplitudes, and £. One 
can impose the on-shell condition 



which implies that the configuration is parity-even. As a result, the field equa- 
tions following from the reduced Lagrangian (2.37) in the static, purely mag- 
netic case are given by 



where N and V(4>) are the same as in (7.13). The boundary conditions are 



w = ±(1 -br 2 ) + 0(r 4 ), = cr + 0(r 3 ), m = 0(r 3 ), (7.35) 
w = -l + B/i exp(-r), = 1 + Cf 2 exp(-er), m = M + O^ 1 ) (7.36) 



at the origin and at infinity, respectively, where the notation is the same as in 



Let us mention the following important differences with the triplet EYMH 
case. First, (7.36) implies that the solutions are neutral. Secondly, the vacuum 
manifold for the Higgs field is now three-sphere, and since ^(S* 3 ) is trivial, 
no topologically stable solutions exist. However, as was argued in [237,208], 
the non-triviality of ir 3 (S 3 ) can lead to the existence of saddle point solutions 
in the flat space YMH theory, which relate to the top of the potential barrier 
between the topological vacua. The corresponding argument is essentially the 
same as the one presented in Sec. 3. 2 for the BK sphalerons. First, one finds 
sequences of static fields {A [A], $[A]} interpolating between the distinct vacua. 
Identifying the end points, they become non-contractible loops. Second, for 
each loop one computes the maximal value of energy, E max , which is then 
minimized over all loops to obtain E s = inf{_E max }. Finally, it is assumed 
that a) E s exists and positive; b) there is a loop whose E m3iX is equal to E s . 



(7.30) 




(7.15). 
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Fig. 13. On the left: Amplitudes w and 4> for the fundamental EYMH sphaleron 
solution (solid lines) and for its first excitation (dashed lines) for k = 0.04 and 
e = 0. On the right: the ADM mass M and the parameter b in (7.35) versus k 
for the fundamental and the first excitation branches. The solid and dashed curves 
merge at the critical value K max = 0.281. The curves for c(k), where c is the third 
parameter of the solutions in (7.35), merge in a similar way. 

Even though these assumptions have not been rigorously proven, they seem 
to be plausible in the case under consideration. This implies the existence of 
a saddle point solution called sphaleron. 

The sphaleron solution to (7.31)-(7.34) in the flat space limit was found nu- 
merically in [103,46]. This corresponds to the choice of the plus sign in (7.35) 
and the gauge amplitude w interpolating between the values one and minus 
one - very much similar to the n = 1 BK solution. The Higgs field amplitude 
is a monotone function. The existence of this solution was established in [76]. 
For small e there is only one solution for a given value of e. For large e there 
are additional solutions for which the amplitude £ is non-trivial, these are 
called deformed sphalerons [212,338]. Sphalerons are unstable. The instability 
is associated with the dynamics of the field £. 

Although unstable, flat space YMH sphalerons admit gravitating generaliza- 
tions [149]. First, for small values of k there is a fundamental branch of regular 
solutions that admit the limit k — > 0. These are qualitatively similar to the 
flat spacetime sphaleron. Next, there is a branch of excited solutions with 
large mass, which also look similar. In particular, the amplitude w still has 
one node (see Fig. 13). However, the excitations do not admit the flat space 
limit and for k — > reduce to the rescaled n = 1 BK solution - exactly in the 
same manner as for monopoles. The excited solutions can be thought of as the 
BK sphalerons sitting inside the YMH sphalerons. Finally, there are solutions 
with higher node numbers. For even values of n one chooses the minus sign 
in (7.35). Similar to the n — 1 solutions, those for n > 1 also exist in pairs, 
each pair containing two solutions with the same nodal structure but with 
different masses. For k — > the n-th solution with lower mass reduces to the 
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n-th rescaled BK solution, while the one with larger mass - to the n + 1-th 
rescaled BK solution. The change of the node number then arises as follows: 
for k — > the field configurations shrink and the amplitudes w with different 
node numbers become practically identical. For example, the amplitude w for 
the n — 1 excited solution shown in Fig. 13 coincides for k — > with the one 
for the n = 2 solution, for which one has w(0) = —1. 

For all n the solutions exist only for a finite range of the parameter < k < 
/t max (n, e). The limiting behaviour for k — * /t max is different from that for 
monopoles. Notice that the RN solution (7.25) does not fulfill the equations in 
the doublet case. As a result, the solutions cannot have a throat connecting an 
interior region and the exterior RN part. What happens instead is that the two 
solutions with the same n coalesce for k = K max , at which point all their pa- 
rameters coincide (see Fig. 13). For example, the two solutions shown in Fig. 13 
get closer and closer as k increases and finally merge for « max (l,0) = 0.281 
becoming one solution. The critical values K max thus correspond to bifurcation 
points. It is worth noting that, unlike the situation in the monopole case, the 
limiting solutions for all n are perfectly smooth and regular. In particular, the 
minimal value of N is quite large. No regular solutions exist for k > « max - 

For any n the regular EYMH sphalerons can be generalized to include a small 
black hole inside. This is possible for < k < K, max (n,e) and for the event 
horizon radius being bounded from above by some value fh(n), which tends 
to zero for k — > K max . For a given n and r h < r h (K) there are again two different 
solutions, which coalesce as r h — > f ft (/t). For r h > r h (n) only the Schwarzschild 
solution is possible. All known regular and black hole solutions of the doublet 
EYMH theory are unstable [47,334,241]. 

The picture described above remains qualitatively the same for any e < oo. For 
e — > oo the Higgs field becomes frozen, = 1, and, if the phase of the Higgs 
field is still fixed according to (7.30), one obtains a simpler theory of gravita- 
ting massive non-Abelian gauge field. This is sometimes called the Einstein- 
non-Abelian-Proca model. [149,231,308,301]. The conical singularity at the 
origin for regular solutions can be avoided in this case by choosing the minus 
sign in (7.35). As a result, solutions for even values of the node number n are 
globally regular. For odd values of n solutions contain the conical singularity, 
as in the monopole case. All regular solutions for r > 0, as well as their black 
hole counterparts, look qualitatively similar to those for e < oo. In particular, 
they exhibit the same bifurcation picture. This model has also been studied 
in the context of the Brans-Dicke theory [302]. 

If the constraint (7.30) is abandoned then the Higgs field is allowed to rotate 
in the internal space spanning a unit sphere S 3 . The EYMH theory reduces 
then to the gravitating gauged 0(4) sigma-model, whose solutions are called 
local textures. The equations of motion in the static, spherically symmetric 
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and purely magnetic case can be obtained by varying the reduced Lagrangian 
(2.37) with = 1. Note that for £' 7^ one cannot set to zero the radial 
component of the gauge field a r . We are unaware of any curved space solutions 
in this case, while those without gravity were studied in [212,338]. For global 
gravitating textures static solutions with finite energy cannot exist, however, 
there are non-trivial time-dependent solutions [120], which can be important 
in the theory of fluctuations of cosmic microwave background [310]. 

7.4 Gravitating Skyrmions 

The gravitating Skyrme model provided one of the first indications of the 
existence of hairy black holes. In [229] it was found that the equation for the 
chiral field on the Schwarzschild background admitted a regular solution. This 
can be regarded as approximately describing a black hole with Skyrme hair. 
In [228] the self-consistent gravitating problem was treated numerically, which 
work for some reasons remained almost unknown. In [147] self-gravitating 
Skyrmions with higher winding numbers were studied as candidates for soliton 
stars, but no stable solutions were found. The first systematic investigation of 
the problem was undertaken in [117,168,167,172]. This revealed the existence 
of stable regular gravitating Skyrmions and Skyrme black holes. The solution 
space of the models was studied also in [38]. The thermodynamics of Skyrme 
black holes was considered in [306]. 

The equations of motion of the Einstein-Skyrme model following from the 2D 
Lagrangian (2.40) read 

m' = N + sin 2 X ) X ' 2 + (r 2 + \ sin 2 X ) (7.37) 

a^Kafr + l sin 2 x)x' 2 , (7.38) 

(vN (r 2 + 2 sin 2 x ) X ')' = o (l + N X ' 2 + ^) sin 2 X , (7.39) 

where N = 1 — 2nmjr and k = AnGf 2 . The length scale is given by L = 1/ef, 
and the dimensionful energy is (Anf/ejM with M = m(oo). The function a 
can be eliminated from the equations. The boundary conditions for the regular 
solutions read 

x = 71-1/ - br + 0(r 3 ), 2m = b 2 (l + b 2 ) r 3 + 0(r 5 ), (7.40) 
x = ar~ 2 + 0(r~ 3 ), m = M + 0(r" 3 ), (7.41) 

at the origin and at infinity, respectively. Here b, a, and M are free parameters 
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Fig. 14. The v = 1 fundamental Skyrmion (solid lines) and its excitation (dashed 
lines). On the left: the chiral field x( r )/ 7T an< i the metric function N(r) for the two 
solutions with k = 0.04. For k — > K max = 0.0437 the solid and dashed curves merge. 
On the right: the mass M and the parameter b in Eq.(7.40) versus k for the lower 
and upper branches of solutions. 

and v is integer. 

The flat space Skyrme model can be regarded as an effective theory in the 
low energy limit of QCD [285,286,336]. The particle-like solutions are inter- 
preted as baryons, and the integer parameter v in (7.40) then plays the role 
of the baryon number. Note that v has the meaning of the topological wind- 
ing number. This is because the description of the chiral field in terms of the 
SU(2) valued function U leads to the mapping S* 3 — > SU(2), provided that the 
asymptotic condition (3.11) is imposed. The winding number is then given by 
Eq.(3.12), v = k[U], which formula reduces in the spherically symmetric case 
to v — {x(0)— x(oo)}/7r. For v — 1 the flat space solution was obtained numer- 
ically in [180]. This solution is stable. Solutions for higher winding numbers are 
also known, but they turn out to be unstable. There are also stable solutions 
for v > 1, but these are not spherically symmetric. It is worth mentioning 
the role of the higher order term F 2 in the Skyrme Lagrangin in (2.38). Due 
to this term, the theory contains both attractive and repulsive interactions, 
which leads to the existence of non-trivial solutions. 

Consider solutions to (7.37)-(7.39) for k ^ 0. The numerical analysis reveals 
that the gravitating Skyrmions exhibit the same typical features as the so- 
lutions of the EYMH models [117,172,38,306]. In particular, since Skyrmions 
have no charge, their behaviour is very similar to that for the EYMH sphale- 
rons. For small values of k there is the fundamental branch of solutions with 
v — 1 that reduce to the flat space Skyrmion as k — > [117]. For these solu- 
tions x, m i an d o are monotone functions, while N develops one minimum at 
some r m . These solutions are stable [168]. Let us call this branch of solutions 
lower branch. Next, for v — 1 there are also excited solutions with larger mass, 
which look qualitatively similar (see Fig. 14) [38,172]. In particular, x is still 
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monotone. For k — > the excitations become infinitely heavy and, remarkably, 
similar to the situation in the EYMH case, reduce to the rescaled BK solution 
with n = 1. Specifically, under 

r —>■ y/nr, m —>■ m/^/K, cosx~^w (7.42) 

the system (7.37)-(7.39) reduces in the limit k — > to the EYM field equa- 
tions, whose solutions are the BK sphalerons. Let us call the branch of excited 
solutions upper branch. 

The lower and upper branches for v — 1 exist only for k < K max . In the 
limit k — > K max the two branches coalesce, exactly in the same way as this 
happens for the EYMH sphalerons. No regular solutions exist for k > K max . 
At the bifurcation point the stability behaviour of the solutions changes: all 
solutions of the upper branch are unstable. This change of stability agrees with 
the general considerations based on the catastrophe theory [306,231,308,301]. 
For higher winding numbers, v > 1, the picture is essentially the same: there 
are two branches of solutions merging at some K max (z/). The value ft m ax(^) 
decreases rapidly with growing v. For k — > the lower branch reduces to the 
flat space Skyrmion with v > 1, while the upper branch again tends to the 
rescaled BK solution; more details can be found in [38,172]. All solutions with 
v > 1 are unstable. It is unknown whether the lower branch solutions with 
v > 1 admit stable and non-spherically symmetric generalizations. 

All solutions described above have black hole analogues. These exist for k < 
K max (y) and if the event horizon radius does not exceed a maximal value 
ffc(/c). The lower and upper branches of solutions exist in the black hole case 
too, and reduce to those of the regular case for — > 0. In the opposite 
limit, rh — > fh^K), the two branches coalesce. No black holes apart from the 
Schwarzschild solution exist for > r^K). The maximal value f/^/c) reduces 
to zero as k tends to K m a X (^)- 

To recapitulate, the critical values K max (z/) for the regular solutions and fh{^) 
with k < K max (z/) for black holes correspond to bifurcation points where the 
lower and upper branches merge. These values determine the region of the 
parameter space («, r^) for which non-trivial solutions exist. A similar bifur- 
cation picture arises in the doublet EYMH theory. It is worth noting that the 
solutions specified by the critical values are perfectly regular. In particular, 
the minimal value of N is non-zero. In the monopole case, on the other hand, 
the limiting solutions exhibit special geometrical structures^]. 

Note finally that, apart from the solutions described above, for which \ is 
monotone, there are also those for which x oscillates around the value n/2. 

17 The charged EYMH black holes with k < 3/4 show the usual bifurcation picture. 
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These exist for very small values of k [38,172]. The solutions of the gravitating 
SU(2) Skyrme model have been generalized to the gauge group SU(3), in which 
case their essential features remain the same [203]. The critical collapse of 
Skyrmions was considered in [39]. 



8 Concluding remarks 

We have given a fairly complete account of known up-to-date solitons and 
black holes in the basic gravity-coupled models with gauge and/or scalar fields 
which respect the non-Abelian symmetries. These solutions show that some 
of the fundamental concepts of the Abelian Einstein-Maxwell theory do not 
extend to the non-Abelian case. Instead, a number of new typical features 
arise. One of such features is the appearance of the characteristic discrete 
structures due to the interaction of the YM field with gravity. These manifest 
in the BK solutions and their various generalizations. For those solitons which 
exist in the flat spacetime limit, one observes the spectrum of gravitational 
BK-type excitations. Another typical phenomenon, which arises in theories 
containing a length scale other than the Planck scale, is the existence of the 
upper bound for the gravitational coupling constant k and for the event hori- 
zon radius rv Non-trivial solutions can exist only for small values of k and 
r/j. For those approaching the bound, the solutions either become gravitation- 
ally closed or develop bifurcations, where different branches of solutions merge. 
Such a merging is accompanied by the change of the stability properties, which 
can be interpreted in terms of catastrophe theory. 

The existence of static and non-spherically symmetric solutions, so far demon- 
strated at the non-linear level only in the EYM theory, is probably also typi- 
cal for gravitating non-Abelian models. This is supported by the perturbative 
analysis for the quite general non-linear model (7.28) carried out in [273]. 
Such deformed solutions can be regarded as the limit of multi-soliton/black 
hole configurations when the separation tends to zero. 

A generic feature associated with the the interior structure of non-Abelian 
black holes is the absence of Cauchy horizons. Let us remind in this connection 
the basic idea behind the mass-inflation scenario [262]. The Cauchy horizon 
inside a Kerr-Newman black hole, unlike the event horizon, exists solely due to 
the high symmetry of the solution. In the generic situation, when deviations 
from the strict axial symmetry are allowed, inner horizons do not exist. Fol- 
lowing this idea, one can argue in the same spirit that an equally efficient way 
to get rid of Cauchy horizons is to add into the system more general matter. 
This is confirmed by the examples of the non-Abelian black holes. 

For the sake on completeness, let us mention briefly also some other important 
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solutions for gravitating Yang-Mills fields which remained outside the main 
text, because they are not of soliton or black hole type. Most of these are 
obtained in the pure EYM theory. 

8.1.1 EYM cosmologies 

The basic idea is as follows. For the SO(4)-invariant spacetime metric 

ds 2 = a 2 (t){dt 2 -dr 2 -sin 2 r(d$ 2 + sin 2 $d<p 2 )}, (8.1) 

with r e [0, 7r], one can find the SO(4)-invariant YM field. This can be obtained 
either as a solution of the symmetry conditions (2.10) [160] or by requiring 
that the stress tensor for the spherically symmetric gauge field (2.13) respect 
the SO (4) symmetry [135]. The result is 

A = i 1 XJdXJ- 1 , where U = exp(irT r ). (8.2) 

The YM equations reduce to 

w 2 + (w 2 - l) 2 = S, (8.3) 

where £ is an integration constant. In view of the conformal invariance, the 
conformal factor of the metric, a(t), does not enter the YM equations. For 
any solution to (8.3) the stress tensor has the perfect fluid structure with p = 
3p = 3£/2a 4 , which allows one [80,160,174,284,283,135] to solve the Einstein 
equations. 

The above solutions can be generalized to spatially open and spatially flat cos- 
mological models [160,135], as well as to higher gauge groups [31,250,251,276]. 
The corresponding superspace quantization was considered, for example, in 
[30,79]. Since the metric in (8.1) is conformally flat and the YM equations 
are conformally invariant, Eq.(8.3) gives rise also to the flat space solutions 
for the YM field describing collapsing shells of the non-linear YM radiation 
[230,122,145]. The homogeneous and anisotropic EYM cosmological models 
were considered in [101,102,16], in which case the solutions turn out to be 
chaotic. The inhomogeneous and spherically symmetric case was considered 
in [282]. 

8.1.2 Cosmological sphaleron 

Eq.(8.3) describes a particle in the double-well potential V(w) = (w 2 — l) 2 . 
For £ — 1 there is a non-trivial static solution for the particle sitting on 
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the top of the barrier separating the two wells, w = [174]. This solution was 
considered in [144,108] as a model of collapsing BK solitons. Since the bottoms 
of the well, w — ±1, correspond to pure gauge fields with different winding 
numbers, which is obvious from (8.2), the static configuration with w = 
can be naturally interpreted as sphaleron [144,107]. This sphaleron solution 
is distinguished by the fact that it consists of the pure gauge field alone, 
which is possible due to the interaction with the background gravitational field. 
Another characteristic feature of the solution is its high (SO (4)) symmetry. 
This implies that the configuration is not localized and fills the whole universe. 

Such a sphaleron solution can be used to evaluate the fermion production 
rate in a closed universe. This problem was considered in [320,322] in the 
approximation when the universe described by the metric (8.1) is static. It 
was assumed that the universe is filled with (quazi)-thermal YM quanta cor- 
responding to excitations over the trivial YM vacuum with w — 1, while the 
sector with unit winding number is empty. As a result, there is a net diffu- 
sion of the field modes between the two sectors, which is accompanied by the 
change in the fermion number due to the axial anomaly. The corresponding 
diffusion rate was computed in [320,322] at the one-loop level. 



8.1.3 EYM instantons 

When continued to the imaginary time, the cosmological solutions described 
above fulfill the Euclidean EYM equations [313,175,314,31,341,270,110]. The 
metric (8.1) then becomes conformal to the standard metric on S 4 , while the 
potential V(w) in the YM equation (8.3) changes sign. Solutions with finite 
Euclidean action are interpreted as tunneling geometries leading to the cre- 
ation of baby universes. Note that, unless the constant S in Eq.(8.3) vanishes, 
the YM field is non-self dual. In the self-dual case one has F^ v = implying 
that the stress tensor is zero. As a result, in order to obtain EYM instantons 
with self-dual gauge fields one can start from a vacuum gravitational instanton 
and use it as a fixed background on which the YM self-duality equations are 
solved. Since topology of gravitational instantons can be quite arbitrary, so- 
lutions for self-dual YM fields can differ substantially from those in flat space 
[84,85,264,119,48-52,81]. It was shown in [337] that the self-duality equations 
for the YM field in axially symmetric case in flat spacetime are equivalent to 
the Ernst equations in GR. 

The results mentioned above, together with those discussed in the main text, 
almost exhaust the list of 4D solutions for gravitating YM fields. Relatively 
little is known about solutions with symmetries other than spherical. At the 
same time, we do not discuss here solutions for D^4, since a separate review 
would be necessary for this (see [106,216,65] for some results in D = 3). A large 
number of such solutions can be obtained from the ten-dimensional heterotic 
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five-brane via various compactifications [118]. 
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